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Preface 


A large amount of ihe recent technical information in the field of forced 
flow and once-through boiling heat transfer and pressure drop has resulted 
from the NASA-AEC Rankine cycle space power programs. This infor- 
mation is in the form of experimental, semi-empirical, and analytical 
functions. A range of working fluids including liquid metals, organics, 
and cryogenics in addition to water and refrigerants has been studied. 

As a result of this work, a number of advances have been made in forced 
flow boiler design. The use of boiler tube inserts to improve heat transfer 
performance is of particular note. 

A number of new forced flow boiler applications have become important 
in the last several years (for example, in mobile propulsion systems, 
topping cycles for power plants, bottoming cycles for power plants, ad- 
vanced refrigeration systems, boiling nuclear reactors, and geothermal 
power systems). There are also many familiar forced flow boiler appli- 
cations where further design optimization may be of value (for example, 
in auxiliary power plants, solar evaporation systems, and process heat 
exchangers involving change of phase). Therefore, this monograph has 
been written with the thought that the newly- developed boiler technology 
may be of value to such boiler applications. 

Experimental, semi-empirical, and analytical boiling information has 
been compiled in this monograph, together with boiling regime transition 
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criteria. These boiler design functions are then applied in example 
boiler designs; the illustrations give the reader sufficient information 
to demonstrate the value of using boiler tube inserts and developing 
optimum boiler designs. 

The monograph contains over 100 technical references in the boiling 
field (many of which contain other references and compilations of other 
contributions which can be found in the literature). It is clear that 
reference omissions exist because of the breadth of the field and space 
limitations. The references presented are thought to be representative 
and useful to those who may want further details. 


Solana Beach, California 
December 1973 


Heinz F. Poppendiek 
Cullen M. Sabin 
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CHAPTER 1: INTRODUCTION 


This monograph contains representative boiling heat transfer and pressure 
drop information obtained primarily from past NASA and AEC programs; 
the results presented are directly applicable to forced flow and once- 
through boiler systems found in many technological areas. 

The primary purpose of this monograph is to present representative boil- 
ing data that desc’ ibe the technology of forced flow and onee-through boiling. 
The information given in this document is intended to be useful for both 
technical and non-technical readers. Chapters 1, 2, and 6 have been 
primarily prepared for managers and those who are not directly involved 
in technical activities. Chapters 3, 4, 5, and the Appendix have been pre- 
pared for technical people such as designers of power systems. 

The forced convection boiler has a number of advantages that should be of 
particular interest to designers in modern technology. In a properly designed 
system, there is little possibility of flow mal-distribution such as sometimes 
occurs in natural circulation systems. Heat transfer characteristics are 
usually consistent, and conductances are predictable, so that higher heat 
fluxes may be employed with safety. This leads to more compact, lighter 
weight equipment. It should be noted, however, that in gas-fired systems 
particularly, the controlling heat transfer resistance may be on the hot side, 
so that increased fluxes would require extended surfaces. If in a power 
generation system the working fluid is very expensive, a forced flow boiler 
can be designed especially for small holdup volume. If the fluid is tempera- 
ture sensitive, the boiling side wall temperatures can be "tailored" to 
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maintain maximum heat transfer rates without overheating the fluid. 
Finally, the forced flow and once-through configurations may be the only 
type which can satisfy a specific need (such as the automotive Rankine 
cycle power plant design having a very short time-response boiler) . 

Forced convection vaporization processes to which the material in this 
monograph is appropriate occur in many areas of technology; in power 
generation, in refrigeration machinery, in vaporization of cryogenic fluids, 
and many other applications . However, since the majority of the infor- 
mation presentad here was originally obtained in Rankine cycle power 
system development programs, the terminology appropriate to power 
systems engineering has been used throughout most of the text. 
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BOILER TYPES 


The simplest form of boiler is a pressure vessel containing a vaporizing 
fluid, with a heat source surrounding it. Usually, however, the heat 
transfer area is increased by utilizing arrays of tubing together with 
one or more larger pressure vessels. 

There are two general classifications of boilers with tube matrices to 
increase heat transfer area. These are the "fire-tube" type, in which 
the vaporizing fluid surrounds tubes through which a heated stream flows; 
and the "water-tube" type, in which the heated stream flows over the 
outside of the tubing and the vaporizing fluid flows within. In each of 
these two classifications, there are also two subclassifications: the 
vaporizing fluid may be circulated by gravity through density difference 
("natural convection") or the vaporizing liquid may be circulated by a 
pump ("forced convection") The distinction between the latter two cate- 
gories is clear from a systems standpoint; however, the boiling 
mechanisms are probably little affected by the means of circulation. 

The once-through boiler is a specific type of forced convection "water 
tube" boiler. 
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HISTORICAL INFORMATION 


It is thought appropriate to present to the reader some historical infor- 
mation on boilers that have been used in technology. This summary is 
not intended to be all-inclusive, but only illustrative of some of the con- 
tributions made in the past. The classifications into which the boilers 
described in the following paragraphs fall are outlined below. 

A. Vaporizing Fluid Flowing Over Tubes ("Fire Tube") 

1. Natural Convection 

a. Horizontal return-tubular 

b. Vertical tube internally fired 

c. Scotch Marine 

2. Forced Convection 

Large air conditioning water chiller evaporator 

B. Vaporizing Fluid Flowing Inside Tubes ("Water Tube") 

1 . Natural Convection 

a. Longitudinal drum 

b. Two drum 

2. Forced Convection 

a. Recirculating types 

(1) Steam pumping (Loeffler) 

(2) Water pumping (LaMont) 
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b. Once through 

(1) Flash boiler (Serpollet) 

(2) Steam automobile (White) 

(3) Mercury vapor Rankine cycle (SNAP-8) 

Horizontal Return-Tubular Boiler 

Figure 1-1 shows the primary elements of a horizontal return tube boiler 
wherein hot combustion gases flow through tubes to heat the surrounding 
water. In this particular boiler, the combustion chamber is external to 
the boiler itself. This type has been used for small, stationary steam 
plants . 

Internally- Fired Vertical Tube Boilers 

Figure 1-2 shows a typical internally- fired vertical tube boiler. In this 
system, the combustion chamber is directly below the vertical fire tube 
array. The type is referred to as internally fired because the lower 
portion of the water tank surrounds the furnace. It is also used for 
small, stationary steam plants . 
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Figure' 1-2. Vertical tube internally filed boiler 




Scotch Marine Boiler 


Geometrically, this boiler is not much different from the horizontal 
return-tub Jar type, except that it is internally fired. The combustion 
chamber is directly below the horizontal boiler tubes, and completely 
inside the boiler shell, with a hot gas passage leading to one end of the 
tube array. These boilers coi tain large quantities of water and steam 
and are utilized in marine ser\ Ice where space is of importance and 
there are advantages to significant energy storage in the boiler. 

Water Chiller Evaporator 

An example of a forced convection boiler in which the vaporizing fluid 
flows over the outside of the tube bank is fou.id outside of the field of 
power generation. The evaporators in many mechanical refrigeration 
machines using cold water as an intermediate heat exchange fluid employ 
such an arrangement. In this case, the water to be chilled flows through 
horizontal tubing instead of combustion gases, and the vaporizing refrig- 
erant is pumped into the tube bank from below . A considerable amount 
of recirculation due to gravity and momentum exchange from the re- 
frigerant inlet stream is also found in a typical evaporator of this type 
since the inlet flow is not usually well distributed along the bank. A 
schematic diagram of such a boiler is shown in Figure 1-3. Not shown 
are the baffles which usually are included to improve flow distribution. 
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Figure 1-3. Water chiller evaporator 




Longitudinal-Drum Water Tube Boiler 

Figure 1-4 shows a schematic diagram of a longitudinal-drum water tube 
boiler. Combustion gases are ducted across the tube bundle by means 
of a baffle system. Typical installations also include a superheater tube 
bank and a water preheater, or n economizer, " which are not shown in 
the diagram. In this boiler, water circulation is insured by the inclined 
tube bank. Direction of the circulation is indicated by arrows. 

Multiple Drum Water Tube Boiler 

There are man}' arrangements for straight and bent tube, two and three 
(or more) drum boilers. Most of these installations include air preheaters, 
water preheaters (economizers) and superheaters. Usually the furnace 
walls are lined with water tubing in addition. 

Steam Pumping Boiler 

In a forced convection recirculating boiler, the designer can choose to 
pump either saturated steam or saturated water. In the Loeffler boiler, 
steam is pumped through a superheater coil, and most of the steam is 
recirculated to the unfired liquid reservoir to produce more steam. This 
system must be operated at high pressure in order to reduce the steam 
volume to be pumped to manageable proportions. The salient feature of 
this boilei is that the liquid supplied to it need not be particularly pure. 
Impurities are retained within the liquid reservoir. A diagram of this 
arrangement is shown in Figure 1-5. 
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Figure 1-5 V:ipor pumping boiler . 


Water Pumping Boiler 

In the LaMont boiler, the designer has chosen to circulate saturate J 
v ater . As in the Loeffler boiler, the liquid drum is unheated, bi t the 
tube surfaces in this case have water pumped through them at se\ eral 
times the rate required to provide the necessary steam flow. Steam 
separation is accomplished in the drum, and the steam produced is 
routed to a superheater. 

Flash Boiler 

A type of once through boiler is the so-called flash boiler, in which the 
preheated vaporizing fluid is sprayed into the hot tubing, partially vapor- 
izing in the nozzle. The liquid droplets film boil under normal 
circumstances, since without elaborate temperature control, the tube 
wall temperature will always be above the transition temperature - 
Some modern work associated with space auxiliary power programs has 
come very near to this design concept. 

Steam Automobile Boiler 

Most steam automobile boilers have been of the once through, single 
coiled tube design. An example of this type is the White Automobile, 
developed in 1907, shown schematically in Figure 1-6. This unit is quite 
similar to some of the boilers being used in steam car developments 
today. One notable difference is that in the White boiler the liquid flow 
was downward through the coil and vertical tubing interconnections were 
required to prevent water from flowing through the coil unvaporized. 
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Figure 1-6. Section through a portion of the White 
automobile boiler. 
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In most of the present day automobile boilers, the first few rows 
of tubes next to the combustor are usually vaporizers, and rows 
farthest from the combustor are devoted to liquid preheating, with 
superheater rows between. 

Mercury Boiler 

An example of a liquid heated boiler is the Space Nuclear Auxiliary 
Power (SNAP-8) mercury system, which utilizes NaK as the heat 
source. The arrangement is pure counterflow, with the vaporization 
taking place in seven identical parallel passages. A schematic draw- 
ing of this system is shown in Figure 1-7. A complete description 
of this unit ' 3 given in Chapter 5. 
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Figure 1-7. NASA SNAP-8 boiler arrangement. 


WORKING FLUID TYPES 


There are a number of working fluids that are used in boiler systems; 
the specific application involved usually defines the fluid type to be con- 
sidered. Liquid metals are normally used in compact, high-power 
density systems because these fluids have excellent heat transfer proper- 
ties and are characterized by relatively low vapor pressure. The latter 
feature means that heavy pressure vessels are not required. Organics 
are also used in power generation systems where low vapor pressures 
are desired in contrast to the high piessure water systems; however, 
the heat transfer characteristics are not as good as water or the liquid 
metals. Freons are used in refrigeration nnd airconditioning systems 
(the evaporators in these systems are special cases of boilers) . Cryogenic 
fluids are normally used in aerospace systems where auxiliary power 
systems are driven by cryogenic vapor bleedoff . 

In Appendix A there are presented in graph form the thermal conductivity, 
specific heat, viscosity, and density of the liquid and vapor phases of 
some of the more practical working fluids found in technology. Vapor 
pressure and surface tension information is also given. 
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CHAPTER 2 


PHYSICAL DESCRIPTION OF BOILING HEAT TRANSFER PROCESSES 


When compared to single phase heat transfer and pressure drop in forced 
flow heat exchange systems, two phase or boiling heat transfer and pressure 
drop are more complex. In boiling flow, the presence of the two phases 
complicates the flow type because the relative location of the two phases 
in a flow tube both in radial and transverse positions is variable. This 
relative location depends on whether the composite flow is linear or helical 
in character, and on the vapor quality and mass flow rate. In this chapter, 
physical descriptions are given of the more important distributions of the 
liquid and vapor phases that can exist in forced flow boiling. These de- 
scriptions are presented as an introductory basis for the experimental, 
semi-empirical and analytical boiling functions given in Chapters 3 and 4. 

The phase distribution descriptions which follow are arranged approximately 
in the order in which they occur along a boiler tube. Figure 2-1 shows the 
general locations of the various phase distributions in a boiler tube in 
relation to the vapor quality (fraction of vapor) and bulk fluid temperature 
profiles. It is pointed out, however, that not all regimes exist in a given 
boiler tube. The criteria that will assist the designer in defining the phase 
distributions in operation are discussed in Chapter 4. It is also reempha- 
sized that because flow conditions change continually along a boiler tube, 
the processes involved are developing rather than steady state. In addition, 
two phase flow has much larger changes in mass flow rates a.,J velocities 
than are found in single phase flow . 
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Figure 2-1. Typical flow regimes in a bone 
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DESCRIPTION OF THE MAJOR BOILING REGIMES 


\ . Nucleate Boiling or Boiling Inception At A Surface 

As a working fluid flows through a boiler, it first encounters a zone 
or component called the preheater wherein the bulk liquid tempera- 
ture is raised to the boiling point. Near the preheater exit the 
temperatures at the tube wall are somewnat greater than the saturation 
temperature (thermodynamic equilibrium boiling temperature). When 
the liquid superheat temperature at the wall reaches a certain critical 
value (depending upon the type of tube material, the number of nuclei 
at the surface, the amount of dissolved gases present and many other 
parameters), ebullition or boiling begins at the tube surface. Bubbles 
are born, grow, and collapse at nucleation sites covering various parts 
of the surface. This process, which repeats itself over and over again, 
is illustrated by an artist's drawing in Figure 2-2. In nucleate boil- 
ing, heat transfer is improved significantly over forced convection 
where no phase change occurs. The improvement is related primarily 
to two mechanisms: namely, turbulent agitation of the boundary layer 

by bubble growth and collapse, * and the added evaporative heat transfer 
2,3 

within the bubble. 

The amount of liquid superheat observed can range up to several 
hundred degrees Fahrenheit, although in most forced flow systems 
it is usually small. Usually bubble generation starts in the region 
where the bulk is still subcooled. Even beyond this point, net vapor 
generation does not necessarily coincide with the location at which 
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the bulk mixed mean fluid temperature is at the saturation tempera- 
ture. The fluid in the center of the tube can still be subcooled wl i'.e 
a stable layer of bubbles is carried along in the wall layers. There 
is frequently a significant void fraction in the flow by the time the 
equilibrium saturation temperature is reached. 

B. Slug Flow 

As one proceeds down the boiler tube beyond the nucleate boiling 
region (in the preheating zone) , a net quantity of vapor is generated 
because the bulk of the liquid is now saturated and bubble collapse 
no longer occurs . Jonsequently, slugs of vapor and liquid exist in 
that part of the boiler tube, each phase moving through the tube at 
about the same velocity. Depending upon the surface tension and 
wetting characteristics of the liquid and the tube wall surface con- 
ditions, the liquid slugs can wet or not wet the tube wall. In the 
former case, large vapor bubbles flow through that part of the tube. 

In the latter case, the flow is characterized by large drops flowing 
through that region of the tube. These features are shown in Figure 2-3. 

C. Annular Flow 

There is another phase distribution type that can occur in the low 
vapor quality region and it is called annular flow. Usually, a low 
velocity liquid annulus exists contiguous to the boiler tube wall with 
a higher velocity vapor core. The opposite arrangement can also 
exist under some circumstances; namely, a vapor annulus at the 
wall with a liquid core. 
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Figure 2-3. Slug How 
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QUALITY 


These two different versions of annular flow {shown in Figure 2-4) 
are controlled by the liquid wetting characteristics, the wall-fluid 
temperature difference and whether vaporization is in operation 
at the solid-liquid interface or the liquid-vapor interface. The 
fluid flow fields in the two phases are related by the common 
velocities and shear stresses at the vapor-liquid interface. 

D. Rivulet and Droplet Flow 

In the intermediate vapor quality region, slug and annular flow distri- 
butions can change to rivulet and droplet flow. When the liquid phase 
mass fraction becomes small, the large drops {which used to fill 
the total tube cross section in slug now) are converted to droplets 
and rivulets that flow along the boiler tube wall: likewise, the liquid 
film {which used to completely cover the tube wail in annular flow) 
is ruptured and converted to rivulet and droplet flow . The surface 
tension characteristics of the system can cause good contact (wetting) 
between the droplets and rivulets with the tube wall or, alternatively, 
vapor film insulation (nonwetting) can exist. These features are 
shown in Figure 2-5. 

E. Entrained Flow 

In all practical boilers, the mix-'d mean vapor quality at that point 
down the boiler tube where the wall is no longer wet is significantly 
less that 100 percent. The remaining liquid must, therefore, be 
entrained in the vapor flow. The quality at which the annular flow- 
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Figure 2-4. Annular flow. 
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ends can vary for over a wide range, with typical values tor good 
designs with wetting fluids ranging from 60 to 80 percent. Under 
extreme conditions with wetting fluids however, the loss of the 
liquid layer on the wall has been shown to occur at quai’ties as low 
as ten percent. 

Some of the major entrainment mechanisms are related to non- 
wetting of the liquid at the boiler tube wall, roll wave action at the 
surface of liquid annuli, spray from bubble generation in thin liquid 
annuli, and vapor generation under a liquid layer at the wall under 
high temperature difference (film boiling) conditions. Fog flow is 
a special case of entrained flow in which the droplets are very small. 
This case occurs in the high quaiity region. When the temperature 
difference between the wall and saturation values becomes sufficiently 
large in forced flow boiling, rivulet and droplet flow' becomes unstable 
in that tiny droplets of the liquid are spalled or shed from the rivulet 
or droplet (see Figure 2-6). When the drops get small enough they 
can even bounce back and forth between opposite sides of a boiler lube 
wall (see Figure 2-7). This disintegration mechanism plays an 
important role in creating fog droplets which comprise fog flow 
(see Figure 2-8), wherein tiny entrained fog droplets are distributed 
across the diameter of a boiler tube with probably the greatest density 
in the central region because evaporation processes are taking place 
near the wall. Small droplets can also be formed by fragmentation 
from large entrained droplets by shear forces. This flow 
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Figure 2-6. Spalling drops. 
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Figure 2-7, Bouncing droplets. 
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Figure 2-8. Fog flow. 
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distribution type occurs in the high vapor quality region as well as 
in the initial part of the superheated region. 

F. Superneated Vapor Flow 

Superheating of the vapor can occur beyond the boiling region prior 
to delivery of the working fluid to the turbine or expander. Sometimes 
a separate component called the superheater is used in the system 
to achieve this step. As indicated previously, small liquid fog droplets 
can, at times, exist in the core of the superheated vapor flow in a 
nonequtiibrium state. 

G- Flo ‘a 4 eg ion Control by Tube Inserts 

In traditional vapor generators, boilers, or evaporators, phase 
separation has beer accomplished by gravity. A few arrangements 
have included centrifugal separators to remove last traces of liquid 
from the vapor flow. Liquid entrainment has usually been minimized 
by employing very long superheater tubing with complex flow paths. 

With the development of space power systems, however, neither 
gravity nor a long, complex flow path could be employed to prevent 
liquid carry-over into the expander, and boiler inserts began to be 
employed. Early boiler tube inserts were simply turbulence enhancing 
devices, although these have been replaced almost universally with 
flow rotating inserts; some of the pioneering work was done at the 
Oak Ridge National Laboratory in 1957. ’ The purpose of the flow 
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rotator is to centrifuge the more dense liquid droplets onto the 
heat transfer surface. 

There are four principal tube inserts which provide lluio rotation 
which have been used in space power plant vapor generators. They 
are described below: 

1. The helically grooved plug 

This device is a rod that has a number of spiral grooves around 
the outside, which fits inside the boiler tube tightly. When this 
device is inserted into the boiler tube, the spiral grooves form 
a number of parallel, small channels lying against the heat transfc 
area where the grooved plug is installed, usually in the preheat 
and low quality end of the boiler. An illustration of the flow- 
pattern type is shown in Figure 2-9. 

2. The twisted tape divider 

This device is formed from a ribbon which is twisted around 

its long axis. When inserted into the boiler tube, it divides 

the tube across the diameter into two channels which spiral 

around the axis. Because of properties of usual boiler materials, 

the pitch-to-diameter ratio of these dividers is limited to about 

7 

fpisr, unless relatively thick ribbons are used. An illustration 
of the flow pattern type is shown in Figure 2-10. 
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3. The helical vane 


This device is a single screw *hread of rectangular cross section 
which may extend all the way from the wall to the cube center- 
line. This geometry can be formed by wrapping rectangular stock 
around a mandrel, or by splitting washers along a diameter, 
welding them edge to edge, and then stretching the resultant helix 
to the proper pitch. In contrast to the twisted divider, the 
helical vane can be made it, any pitch. Frequently, the vane is 
attached to a rod or tube called a eenterbody, which lies along 
the axis of the boiler tube. This element serves to support the 
vane and can carry instrumentation as well. Figure 2-11 illus- 
trates the flow pattern for this insert. 

4. The wire coil 

This device is a helical wire coil which lies against the wall of 
9 

the boiler tube.' Although it is usually mat,. . round wire, it 

can be made of square wire as well, and in that form it is identical 
to a helical vane of square cross section. It can be made in any 
pitch. The resulting flow patterns are shown >n Figure 2-12. 
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Figure 2-11. Helical vane. 
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PHOTOGRAPHS OF SOME OF THE BOILING REGIMES 


A number of photographs have been compiled which illustrate some of 
the phase distributions that have been discussed above: they are shown 
in Figures 2-13 to 2-18- In some cases, simulated forced flow boiling 
is depicted using water-air mixtures in glass pipes. In other cases, 
transparent heating films and optical ports have been used so that actual 
forced flow boiling phase distributions can be seen. 

In the previous discussion of slug, annular, rivulet and droplet flows, 

it was noted that the effect of wetting on the phase distributions must 

be considered. The degree of wetting of a liquid on a heat transfer surface, 

defined in terms ot a contact angle, can vary significantly w ith temperature 
10 

level See Figure 2-19, for example, where the contact angle of a mercury 
drop on a clean tantalum surface in a mercury vapor environment is shown 
to change markedly w ith temperature beyond 700” F. Actual photographs 
shov i> g two illustrative contact angles measured are shown in Figures 

2-20 and 2-21. 
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Figure 2-15. 
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Linear, wetted-wall slug flow (Freon at low vapor quality). 
Courtesy of General Klectrie Company, Schenectady, New York. 
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Fi ~ un ‘ - ■ 1" ■ Spalling mercury droplet on a heated plate 
Geoscience i.td photograph. 



Figure 2-18. Helical forced flow boiling using a 


w ire coil insert (water-air 
simulation) . 

Courtesy General Electric Company, 
Cincinnati, Ohio 
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Figure 2- 19. Contact angle formed between a mercury droplet 
and a chemically etched tantalum surface versus 
furnace temperature. 
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DISCUSSION 


The material presented in this chapter gives a physical description of 
the various ways in which the liquid and vapor phases could be dis- 
tributed in a boiler lube. When compared to single phase heat transfer 
and pressure drop in forced flow heat exchanger systems, two phase, 
or boiling heat transfer, is more complex. One immediate complication 
of two phase flow is that the relative location of the phases in a flow tube 
is a variable. The phase location depends on whether the composite 
flow is linear or helical in character and on the vapor quality and mass 
flow rate of ihe fluid. Quantitative experimental and analytical de- 
scriptions of these boiling heat transfer regimes are presented in 
Chapters 3 and 4 
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CHAPTER 3: REPRESENTATIVE 
EXPERIMENTAL, SEMI- EMPIRICAL AND 
ANALYTICAL INFORMATION ON 
BOILING HEAT TRANSFER AND 
PRESSURE DROP 


In Chapter 2, physical descriptions of the major phioe distributions 
for boiling heat transfer that are known to exist in forced flow were 
presented. The oresent chapter outlines representative experimental, 
semi-empirical, and analytical information on the frictional pressure 
drop* and heat transfer in tubes with circular cross sections for these 
different phase distributions or flow regimes . In most cases, not all 
three types of information are available for a given regime. Whenever 
oossible, however, comparisons between the experimental, empirical, 
and analytical functions are made. These three types of information 
have a range of applicability. The design engineer must consider each 
information set and, on the basis of the limitations given, determine hou 
he can use it in design. The material on transition phenomena in Chapter 
4 will also assist him in establishing where and when these various boil- 
ing regimes are in operation. 

Sufficient information is presented in Chapters 3 and 4 so that reasonable 
boiler designs can be carried out. Examples of such designs and their 
comparison to actual performance are presented in Chapter 5. 


* 

The total pressure drop is the sum of the frictional and momentum change 
terms. The latter quantity is defined in terms of differences between boil- 
ing fluid cutlet and inlet momentum rates for the cases of separated and 
mixed Pmws. 


51 


FORCED FLOW BOILING HEAT TRANSFER 


AND PRESSURE DROP INFORMATION 


The major forced flow boiling regimes found in a practical boiler 
system as outlined in Chapter 2 are presented in this section. 

* 

/ . Liquid Preheating 

Experimental Information 

Classical information on single phase pressure drop and heat 

transfer in terms of friction factors and Nusselt modul’ are 

well-described in the literature. The experimental viscous 

and turbulent flow friction factors are shown in Figure 3-1 for 

tubes with circular cross sections. 1 The experimental Nusselt/ 

Prandtl modulus versus Reynolds modulus function for turbulent 

flow in tubes with circular cross sections is shown in Figure 

2 

3-2 for ordinary fluids (Prandtl number greater than about 0.5) 

O 

and in Figure 3-3 for liquid metals'* (Prandtl number less than 
about 0.05). 

Semi-Empirical Information 

Also shown in Figure 3-2 is the classical semi-empirical 


♦ 

Liquid preheating has been added to the beginning of the boiling regimes 
and vapor superheating added at the end to complete the overall forced 
flow boiling process description. 
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FRICTION FACTOR 



Figure 3-1. Friction factors for fully-established, laminar and 
turbulent single phase flow (circular cross section). 
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Figure 3-2. Nusselt modulus for fully-established, turbulent single 
phase flow of ordinary fluids (circular cross section). 
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ANAlYTICAl, 



SEE EQUATION SET 3-3 IN APPENDIX ft 






convective heat transfer function which coincides with the mean 

experimental function. A specific semi-empirical correlation 

for liquid metal convective heat transfer was developed for 

4 

turbulent flow in tubes with circular cross section. It is 
shown in Figure 3-3. 

Analytical Information 

The classical viscous flow friction factor for circular cross 
section tubes* is also shown in Figure 3-1; agreement with 
experimental data is excellent. 

A number of analytical expressions for forced flow convective 

transfer have been developed in the literature. Nusselt modulus 

functions for viscous flow in tubes with circular cross sections 

range from 3.65 for uniform wall temperatures to 4.36 for 

5 

uniform wall heat fluxes. The most well-known turbulent 

flow Nusselt modulus function for all Prandtl moduli is based 

on the analogy between heat and momentum transfer and the 

existence of the classical generalized velocity profile defined 

0 

by a laminar sublayer, buffer layer, and a turbulent core. 

The results of this work can be seen in Figure >-A. An 

analogy solution for liquid metals in circular cross section 
7 

tubes is also shown in Figure 3-3. 

B. Martinelli Representations for Two Phase Flow 

A -'mplified description of forced flow boiling in extensive use 
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Heat-momentum transfer analogy for a smooth tube 
(applicable to both ordinary fluids and liquid metals) 





was originally developed by Martinelli. Although both fluids are 
postulated to be in the turbulent state, their distribution in the 
tube is undefined. This gross, nondetailed representation of two 
phase flow has been used successfully to estimate pressure drop 
in boiling systems. Experimental water-air pressure drop data 
were generalized by analysis and the results extended to forced 

g 

flow water boiling. Information on the ratio of two phase pressure 

gradients to the pressure gradients for all-liquid flowing alone is 

shown in Figure 3-5. The local pressure drops are integrated 

along the boiler tube to obtain the total frictional pressure loss. 

Similar procedures have been carried out for estimating the 

g 

pressure losses in forced flow sodium boiling and potassium 

.... 10 

boiling. 

C . Slug Flow 

Experimental Information 

Figure 3-6 shows a graph of some experimental boiling mercury 
pressure drop data as a function of vapor quality under non- 
wetting, linear slug flow conditions. ** The data shown have 
been corrected for the acceleration terms so that only frictional 
pressure drop is involved. The corresponding heat transfer data 
are shown in Figure 3-7. These data were obtained in an electri- 
cally-heated tube system (defined by a nearly uniform wall heat 
flux condition) . 
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Figure 3-7. Comparison of experimental and predicted boiling mercury conductances 
for nonwetting slug flow. 



Analytical Information 

A two phase flow model consisting of alternate nonwetting slugs 
of liquid and vapor, each of which fills the cross section of the 
tube, has been proposed. Although the instantaneous flow at a 
cross section is unsteady, the time average mass flow rates of 
the liquid and vapor phases may be used to define a time average 
vapor quality. This model has been analyzed to determine the 
average frictional pressure drop of the flow. 11 The model was 
evaluated for the case where physical fluid properties are invariant, 
the slug friction factors are equal and invariant, and the vapor 
quality varies linearly with axial distance (uniform wall heat 
addition); the total frictional pressure drop function along a boiler 
tube was determined and is shown in Figure 3-6. A corresponding 
heat transfer analysis 11 (which is outlined in Appendix C where the 
method is illustrated), gives a normalized heat transfer conductance 
versus vapor quality along the boiler tube (see Figure 3-8) . 


An analytical prediction of void fraction in wetted wall slug flow 
12 

has been presented. This expression, which helps to define 

the flow state, gives fairly good agreement for adiabatic systems 

of dissimilar fluids such as air and water. It is shown compared 
13 

to some data for steam and water in Figure 3-9. Some of the 
discrepancy between theory and experiment is probably attributable 
to the fact that the experimental system was an annular test section. 
The inner cylinder of the annulus is one plate of a capacitance probe 
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LOCAL VAPOR QUALITY, x 


Figure 3-8 Dimensionless boiling heat transfer conductances for 
the nonwetting laminar and turbulent slug flow models. 





Figure 3-9. Void fraction versus steam quality for upward 
wetting slug flow. 
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from which the void fraction was determined. In a system with 
heat transfer, the void fraction is not zero at zero vapor quality, 
because 'he void fraction in subcooled boiling (a departure from 
thermodynamic equilibrium) is in typical systems as much as 0.5 
(see, for example. Reference 14) . 

D. Annular Flow 

Experimental Information 

A number of experimenters have studied pressure drop and heat 

transfer in forced flow boiling systems wherein the two phases 

distribute themselves in an annular arrangement (usually the 

liquid phase is at the vail and the vapor in the core) . Some 

examples of such data are shown in Figures 2- tO, 3-11, and 
15 1.6 17 18 

3-12. ’ ’ * The relatively large amount of scatter in 

Figure 3-11 suggests that for the high points, both phases were 
probably in the turbulent state or large surface roughnesses 
existed at the liquid layer surface 

Semi- Empirical Information 

Algebraic correlations for turbulent-turbulent annular upward 

flow of air and water based on some theoretical models have 
19 

been established. The correlations include both void fraction 
and two phase frictional pressure drop. These expressions agree 
well with Martinelli's data, which are based on horizontal flow 
of a number of fluid combinations. Therefore, it is plausible 
that the correlations have application to fluids other than air 


65 



VOID FRACTION 


T 



WvWvPjf 

Figure :i-10. Experimental void fraction in upward annular flow-analytical predictions also shown. 




Figure 3-11 . Experimental two-phase pressure drop in upward 
annular flow — analytical predictions also shown. 
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Figure 3-12. Comparison of superheated annular liquid film model with boiling potassium data. 




and water. The pressure drop correlation (shown in Figure 3-13) 
is specifically restricted to cases in which the liquid entrained in 
the gas phase is less than 20 percent of the total local liquid flow, 
and in which viscous forces in the liquid film are negligible. 

Analytical Information 

A number of annular flow pressure drop analyses for different 

flow regimes in the annulus and core regions have been per- 
20 

formed; typical results for viscous-turbulent flow are shown 
in Figure 3-14. An outline of this solution is presented in 
Appendix C. These results are also compared to experimental 
data in Figures 3-10 and 3-11. 

The heat transfer in an annulus system can be defined by the 
temperature drop across the liquid or vapor annulus contiguous 
to the heat transfer wall . In the former case, heat conduction 
through the liquid film supplies heat to the liquid-vapor inter- 
face where boiling takes place at the saturation temperature. 

In the latter case, convection and radiation through the vapor 
film control the heat of vaporization at the liquid-vapor inter- 
face. Examples of three analytical representations are shown 
in Figure 3-12, where predictions are compared with experi- 
mental liquid boiling information. These models are based on 
idealized, annular, liquid layers contiguous to a wall. The fir t 
model consists of a superheated liquid layer that is heated f» >m 
the tube wall and loses heat from the liquid-vapor interface by 
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DIMENSIONLESS SINGLE-PHASE GAS PRESSURE DROP ,Ap 



DIMENSIONLESS SINGLE-PHASE LIQUID PRESSURE DROP,A P| * 

Figure 3-13. Correlation for two-phase pressure drop in annular flow . 
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Figure 3-14. 


Viscous-turbulent annular flow wit.i smoot i interlace 

<£ t -0.03). 
p 


71 




evaporative heat transfer; heat transmission through the layer 
consists of conduction only. The second model is defined by a 
superheated liquid layer with wall heat addition and a uniform 
heat sink representing volume boiling. The third model is de- 
fined by a superheated liquid layer with wall heat addition and 
a volume heat sink that varies linearly from a maximum value 
at the wall to zero at the liquid-vapor interface. The resulting 

analytical Nusselt moduli for these cases can be shown to be 

21 

D/<5 , 2D/3 , and 3D/6 , respectively- The liquid film 
thickness 6 is obtained from Equation Set 3-14 in Appendix B- 

A closed form algebraic relationship between the void fraction 

and two phase pressure drop has been derived for turbulent- 
19 

turbulent flow. The model on which the analysis is based 

utilizes a one-seventh power velocity profile in the liquid 

phase. It is not necessary to define the properties of the core 

flow, or the nature of the interaction between the core flow and 

the annulus to obtain this relationship. A comparison of the 

analytical result and the experimentally-established relationship 
23 

from Martinelli's data is shown in Figure 3-15. 

Prediction of void fraction from the properties of the two streams 
and an interaction model to be applied at the interface has been 
attempted by several authors. This is a difficult task, since the 
models must represent by time invariant functions some processes 
which vary rapidly with time in an apparently random manner. 
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a. VOID FRACTION 


Figure 3-15. Relationship between void fraction and two-phase 

pressure drop for turbulent-turbulent annular flow'. 
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The apparent liquid-vapor interface, which one may compute from 

void fraction data, is irregular, wavy, and unsteady in time. For 

these reasons, predictions of void fraction from physical models 

have not been very successful. An example of predicted void 

22 

fraction and pressure drop, utilizing a flow model in which both 
streams are characterized by one-seventh power velocity profiles, 
is shown ’n Figure 3-15. 

The annular flow heat transfer model based on conduction through 
the liquid film gives good results for laminar liquid layers or for 
liquid metals, where the thermal conductivity is large compared 
to turbulent transfer. However, for thick liquid layers of ordinary 
fluids, the heat flow through the film would be expected to be con- 
trolled by turbulent convection. This case may be treated by 
replacinr the thermal conductivity in the previous conduction heat 
transfer mouel by a mean effective thermal conductivity con- 
structed with the aid of an appropriate turbulence me del. One 
such mean turbulent conductivity, based on the idealizations 
that the velocity profile in the liquid layer is identical to the one- 
seventh power velocity profile measured for a tube running full 
of liquid and that the turbulent Prandtl modulus is unity, is shown 
in Figure 3-16. The liquid-vapor interface radius, r., in Figure 
3-16, is obtained from the corresponding pn. ssure drop analysis. 
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Figure 3-16. Mean effective thermal conductivity of a turbulent annular 
liquid layer of an ordinary fluid as a function of interface 
radius. 
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E. Rivulet Flow -Wetting 

When a liquid film, which is clinging to the wall of a boiler tube, 
ruptures because of some hydrodynamic disturbance, a wetting 
rivulet results that maintains direct contact with the wall. 

Analytical Information 

25 

An analysis for this phase distribution has been made; it 
accounts for vaporization at the liquid-vapor interface of the 
rivulet and convective heat transfer from the bare portion of 
the tube wall. The resulting normalized boiling conductance 
function is shown in Figure 3-17. The contact fraction, F, for 
the rivulet is defined by the contact angle* and the vapor quality. 

This conductance function predicts the heat transfer in a boiler 
tube from the point where a complete liquid annulus exists (with 
only surface vaporization operating) to the point where the liquid 
rivulet on the wall is all gone (with only vapor convection remaining) . 
The cha..ge in heat transfer between these two end points is 
normally equal to several orders of magnitude. 

F. Entrained Flow 
Experimental Information 

Considerable experimental data are available on annular flows 


* 

Contact angles are generally determined experimentally for the par. icular 
surface, liquid and vapor present. 


76 




with liquid entrainment in the core, for speeP *.■ iystems. A flow 
regime map showing the nature of the bouncily tor entrainment, 
taken from Reference 19, is shown in Figure 3-18. A' V-.<ugh 
this figure is for air-water simulated boiling and cann~ oe used 
to predict entrainment in a general case, the regime dependence 
on velocities of the phases may be expected to be qualitatively 
the same for most ordinary boiling fluids (excepting liquid 
metals). Additional data on entrained flows in simulated boil- 
ing may be found in References 26 and 27. The number of 
physical parameters which can potentially influence the pressure 
drop and phase distribution in an entrained flow is large. For 
this reason, some theoretical or semi-empirical treatment is 
necessary to make useful predictions of conditions n a boiler. 

The following sections outline so ne useful procedUL^s. 

Semi- Empirical Information 

Pressure drop in the so-called annular-dispersed flows, which 
are flows with a continuous liquid film on the wall and a significant 
portion of the cotal local li .mid flow entrained as fine droplets in 
the vapor stream has been treated successfully by a -^odel repre- 
senting the flow as a homogeneous mixture of the tu fluids moving 
at a single velocity. 

In Reference 28 it i: shown that such a model, used with a constant 
friction factor of 0.02, gives an acceptable fit to adiabatic pressure 
drop data for a wide range of fluids. The fluid combinations include 
air-water, air-Heptane, Argon -Ethyl Alcohol, and steam-water at 
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Figure 3-18. Boundaries of the annular flow regime for upward 
flow of water and air (from Fv.Vrence 19). 
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a wide range of pressures. This approximation is applicable at 
entrained flows above 20 percent of the local total liquid flow. 

Theie is a significant pressure drop associated with the entrainment 
process, since the liquid droplets must be accelerated to a velocity 
approaching that of the vapor. In Reference 29, it is shown that a 
good estimate of the acceleration pressure gradient is obtained by 
the difference of product of the local entrained liquid mass flow- 
ana the superf.''i.'l vapor velocity at two stations, divided by the 
distance between the two stations. The superficial vapor velocity 
is the local vapor volume flow divided by the tube 'ross sectional 
area. This approximation was originally proposed in Reference 
30. 

Theoretical Information 

An analytical treatment of annular flow with large entrainment is 
presented ir Reference 31. This analysis consideis the momentum 
and mass transfer components of the shear at the liquid-vapor inter- 
face. It is shown that the momentum term controls in the liquid 
film, and the mass transfer term controls in the vapor core. The 
vapor core is idealized to be a homogeneous mixture, and the liquid 
layer is defined by the universal velocity profile. In both regions, 
the mixing length is postulated to he proportional to the distance 
from the channel wall. A semi-empirical correction is employed 
to account for slip in tne coie. A second semi-empirical correction 
is also employed to account for wave height, and empirical data 
are used to evaluate a universal function for the normalized inter- 
face radius. Th’s function, which was evaluated using data for 
water-argon and alcohol -argon, is shown in Figure 3-19. 
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Figure 3-19- Relationship between pressure gradient and liquid 
film thickness for annular flow with entrainment 
<from Reference 31). 
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A prediction system relating the film thickness and pressure 
drop was developed in the same reference for the additional 
restrictions that the gas density is much smaller than the 
liquid density and the liquid velocity is much smaller than the 
vapor velocity. The prediction procedure is as follows: 

1. Compute pressure drop from any method or correlation. 

2. Obtain liquid film thickness from Figure 3-19 and 
Equations (1) and (2) in Equation Set 3-19. 

3. Compute liquid fi'm flow rate using Equations (3), (4), 
and (5) in Equation Set 3-19. 

4. Compute liquid entrainment from total liquid flow and film 
flow rates. 

A comparison of some experimental data with this prediction method 
is shown in Figure 3-20. 

G- Superheated Vapor Flow 

Semi-Empirical Information 

Superheated vapor flow is identical to ordinary gaseous convection. 
However, since superheater lengths are frequently hydrodynamically 
short and the flow in the boiler is a devekpir , one, so that fully- 
established flow may not be obtained in the superheat region, some 
judgment is required in choosing the proper heat transfer relationship. 
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GAS WEIGHT FRACTION 


Figure 3-20. Comparison of predicted film thickness with 

experimental data for annular flow with entrainment 
(from Reference 31). 
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However, since the vapor flow is almost invariably in the 
Reynolds modulus range associated with the turbulent flow 
regime, the superheater conductance may be expected to tend 
toward that given by relationships for fully-established turbulent 
flow . The fully-established correlation can then be modified by 
a term dependent upon normalized tube l eng th . 

Under some circumstances, partic ularly in uniform heat flux 
systems or where cryogenic fluids are being vaporized, the wall- 
to-bulk temperature differences can be large compared to the 
absolute temperature, so that the thermal properties cf the fluid 
may vary significantly through the wall boundary layer. A number 
of author ’ have proposed methods of compensating for variable 
fluid properties. Usually the correction procedures involve the 
use of U’Md properties based on the bulk temperature, the wall 
temperature, a mean "film" temperature, or a combination of 
these properties. By use of some relationships for ideal gases, 
it is possible to reduce these various correlations to expressions 
based on bulk fluid properties and a correction term. Eight 
. orre.ations (References 32 through 39), based on data for air, 
hydrogen, carbon dioxide, and helium, have been reduced to a 
common form, and the correction terms compared. The 
correction terms for tube length are shown in Figure 3-21: 
the correction terms for temperature are shown in 



u 


Figure 3-21 . Change .n Nusselt modulus with tube length 
acco'\' ; "g to several correlations. 
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Figure 3-22. Also shown for comparison is the well-known 

2 

correlation given by McAdams . 

The friction factor is also affected by temperature dependent 
fluid properties and developing flow. The adiabatic fully-established 
friction factor can be modified by multiplicative correction terms 
in the same manner as the Nusselt modulus relationship. The 
temperature variation of a typical friction factor relationship* 
reduced to this multiplicative form, is compared to experimental 
data in Figure 3-23. Another such friction factor correlation 
with a multiplicative wall-to-bulk temperature ratio correction 
is given in Reference 36. This latter correlation is, however, 
based on the Reynolds modulus constructed with properties 
evaluated at the wall temperature. Typical tube length functions 
are similar to the ones used in heat transfer relationships; their 
variation can be found in fluid dynamics texts. 


* 

The relationship shown in Figure 3-23 was derived from an expression 
taken from Reference 35 by use of a Sutherland model for the vapor viscosity. 
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Figure 3-22. Effect of wall to bulk temperature ratio on 
Nusselt modulus according to several 
correlations. 
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Figure 3-23. Temperature variation of friction factoi fur 
single phase vapor flow. 
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H. Flow Through Tubes with Inserts 

The material on flow through tubes with inserts is divided into 
two parts, namely, (1) a general description of helical flow, and 
(2) a discussion of flow in specific insert systems currently in 
use. 

General Characteristics of Helical Flow. 

1. Single Phase Flow 

Experimental Information 

Available single phase helical flow friction factor data in the 

literature (primarily measured under adiabatic conditions) 

40 

have been converted to a common basis. The mean experi- 
mental curves presented in terms of helical to linear friction 

4 

factor ratios as a function of pitch to diameter ratio for 10 
5 

< Re < 10 are shown in Figure 3-24. Helical-flow forced 
£ 

convection data available in the literature obtained for heating 

40 

water, air, and ethylene glycol were also correlated. 

Mean experimental helical to linear Nusselt modulus ratios 

4 5 

as a function of pitch to diameter ratio for 10 < Re <10 

f 

are shown in Figure 3-25. 

Analytical Information 

The increased heat transfer in single-phase helical-flow systems 
results primarily from increased fluid turbulence and free 
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Figure 3-24. Ratio of helical to linear friction factor versus 
pitch-to-diameter ratio. 
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Figure 3-25. Ratio of helical to linear flow Nusselt modulus versus 
pitch-to-diameter ratio. 
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convection (under conditions of wall heat addition) . Exact 
analytical descriptions of the heat and momentum transfer 
in helical flow systems are difficult to obtain because of the 
complex geometries and fluid flow patterns involved. Flow 
instability in curvilinear flow systems and superposed natural 
nvecticn processes further complicate the transport processes 
in helical flow. It is possible, however, to examine the relative 
importance of some of the transport models which can be used 
to characterize single-phase helical flow svstems. 

A number of fluid dynamics analyses applicable to helical or 

rotational flow can be found in the literature. An interesting 

41 

study has been reported for single-phase flow in a tube with 
a twisted divider. An approximate pressure drop expression 
is derived which accounts for axial and tangential flow and 
vortex mixing. The predicted results are in good agreement 
with several sets of experimental data. 

A simplified expression for the isothermal friction factor in 

the turbulent flow regime can be derived by considering only 
40 

gross flow features. A large fraction of a fluid flowing in 

a tune exists in an anuuiar region near the wall. If the fluid 

flows helically through the tube, the effective flow length has 

been increased by the ratio of the helical path at the tube wall 

2 1 1/2 

per unit of axial distance, 1 + ( 7fD/p) J . Similarly, 
the resultant fluid velocity in (ho annular region near the 
wall has increased approximately by the same ratio. If these 


J2 



quantities are substituted into the pressure drop equation, a 
group of terms results which is defined as the friction factor 
for helical flow; this function is shown in Figure 3-24. 

Several heat transfer analyses applicable to helical or 

rotational flow can be found in the literature. One involves 

42 

heat transfer conductances in curved channels. An ex- 
pression for convective heat transfer in a tube with a twisted 

41 

divider was also described . The effect of free convection 
on the overall heat transfer was not included. 

An expression for the turbulent Nusselt modulus for helical 

flow (excluding free convection) has been derived on the same 

40 

basis given for the simple friction factor described above. 

Free convection heat transfer expressions for idealized 
natural convection cells in a horizontal fluid layer under un- 
stable temperature gradient conditions have also been developed. 
Typical power law Grashof-Prandtl modulus equations result 
that are similar to the experimental expressions for turbulent 
free convection from horizontal cylinders. The sum of such 
analytical normalized forced and free convection conductances 
in terms of Nusselt moduli can be used to estimate single phase 
helical flow; a comparison with experimental information is 
shown in Figure 3-25. 
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2. Helical Vapor and Droplet Flow (Nonwetting) 

Consider the case where droplets and vapor ire flowing 
helically through a boiler tube wherein the droplets flow 
along the wall but do not wet it. The heat transfer from 
the tube wall consists of both gaseous convection and 
droplet vaporization. Because these two mechanisms are 
operating in parallel, the boiling conductance for the tube 
is equal to the sum of the gaseous convection and droplet 
vaporization conductances . 

Analytical Information 

An analysis for the droplet vaporization conductance has 
25 

been made and the results in the form of a normalized con- 
ductance as a function of vapor quality are shown in Figure 
3-26. This model describes film boiling vaporization from 
beneath a droplet contiguous to a heat transfer surface; heat 
conduction through the vapor film under the drop causes drop 
evaporation. The centrifugal force field on the drop is 
balanced by the pressure field under the drop. The vapor con- 
ductance (which is the remaining part of the boiling conductance) 
for helical vapor flow was described previously in this chapter. 

A comparison of predicted and experimental potassium boil- 
ing conductances for this type of phase distribution can be 
seen in Figure 3-27. 
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Figure 3-26. Normalized evaporative conductance as a function of 
vapor quality for nonwetting rivule. flow with helical 
insert. 
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Figure 3-27 . Predicted and experimental helical flow boiling potassium results. 
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Flow Through Specific Inserts. 


(1) Helical Grooves in a Rod 

If helical grooves are machined in the periphery of a rod 
that is inserted in a boiler tube, boiling heat transfer 
takes place in an efficient manner. The two phases are 
readily separated by the centrifugal force field with a 
thin liquid layer located at the boiler tube wall in the 

43 

grooves. An analysis of this flow system has been made 
and the details are found in the Appendix. 

(2) Twisted Tape Insert 

When the fluid flowing through a straight tube is brought 

into rotation by a twisted divider across its diameter, 

the flow is forced into a path and configuration that has 

an equivalent diameter, cross sectional area, and velocity 

different from those in the original tube. A mathematical 
44 

relationship for the equivalent diameter of this insert 
geometry is shown in Figure 3-28. 

(3) Helical Vane 

If the twisted divider discussed in the previous paragraph 

only extends halfway across the tube diameter, a helical 

vane system results. The corresponding equivalent diameter 
44 

expression for this geometry is also shown in Figure 3-28. 
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?jp,’:re 3-28. Equivalent diameters for tubes with helical 
inserts (from Reference 44). 
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Helical Wire Coils. 


The preceding sections contain a discussion of heat transfer and 
pressure drop in tubes in which the flow follows a helical path. 

In tubes containing twisted dividers, grooved plugs, or helical 
vanes which extend nearly to the tube centerline, the helical path 
is well defined. This is not necessarily the case for helical wire 
coils. There is, of course, no substantial difference between a 
coil made up of wire equal to the tube radius and a helical vane the 
height of which is equal to the tube radius . However, in a number 
of cases wire coil inserts have been used in which the wire diameter 
is small compared to the tube radius, so that the bulk of the fluid 
can slip with respect to the wire. The helical path of the fluid can, 
therefore, be significantly different from the path defined by the 
wire helix. 

Wire coils can be made up from wire of any cross section. However, 
in the following sections, the data and analyses pertain specifically 
to coils made from round wire. 

Experimental Information 

Experimental programs to determine the heat transfer and pressure 
drop characteristics of single phase flow through tubes containing 
wire coil inserts are reported in References 37 and 45. Results 
of those investigations are shown in Figures 3-29 and 3-30. 

Semi- Empirical Information 

For very short pitch helical wire inserts, the slip of the flow with 
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Figure 3-30. Ratio of measured pressure drop for wire coils to 
that calculated for axial flow in a smooth tube 
(adapted from Reference 45). 
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respect to the wire rotation would be expected to be large, so that 
the flow is similar to that in a tube containing annular rings, or 
corrugations. The limiting minimum pitch for a helical wire coil 
is that pitch at which adjacent coils touch each other. At this 
condition the interior wall of the flow passage appears as a corru- 
gated wall tube whose minimum diameter is the inner diameter of 
the wire coil. 

In Figure 3-30 it was seen that the friction factor passes through a 
maximum. This maximum apparently occurs at the pitch at which 
flow reattachment just occurs between the wires. Very little 
flow rotation would be expected under these conditions . A semi- 
empirical correlation of the pressure drop characteristics of short 
pitch helical wire inserts is presented in Reference 25. The 
results of this work are compared to experimental data for both 
helical wire coil inserts and corrugated tubes in Figure 3-31. 

In the analysis of Reference 25 , it is shown that for short pitch 
wire coils in which the flow reattaches between coil turns, the 
dominant term in the pressure drop equation is the drag of the wire 
coils so that the friction factor should vary inversely as the wire 
coil pitch; i.e., the friction is directly proportional to the number 
of drag elements per unit of length. That this is the case may be 
seen in the figure. 

Analytical Information 

An analysis of the friction and flow rotating characteristics of long 
pitch (p/D >1), wire coils in fully-established single phase flow is 
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presented in Reference 25 . The mathematical model is based upon 
solid body rotation for the fluid, a one-seventh power velocity 
profile, and wire drag and wall friction coefficients which are 
independent of the Reynolds modulus . The pressure drop results 
of this analysis are compared to experimental data in Figure 3- 
32 • The predicted rotation of the flow is shown in Figure 3-33. 

In this same reference, the tube length required for the flow to 
attain equilibrium rotation is also considered. This length is 
such that for practical boilers and wire coil pitches fully- 
established flow will probably not be attained. This topic is con- 
sidered in more detail in Chapter 4, which takes up flow 
transitions . 
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WITHOUT SUP 



Figure 3-33. Rotating component of velocity as a function of 
wire coil pitch and wire diameter. 
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DISCUSSION 


This chapter has presented representative quantitative information on 
the important phase distributions that define forced flow boiling heat 
transfer and frictional pressure drop. The parameters vapor quality, 
liquid wetting, fluid flow rate and insert type, among others, control 
the phase distributions. The qualitative boiling conductance functions 
as a function of vapor quality are shown in a schematic graph in 
Figure 3-34 . The wetted wall annular flow conductances are the 
greatest (because of the existence of very thin vaporizing liquid film 
at the wall) and the nonwetting slug flow conductances are the smallest 
(because of the controlling, poor vapor convection) . The conductances 
for these two flow regimes can, at times, differ by three or four 
orders of magnitude. 

In the higher vapor quality regions, linear flow yields low boiling con- 
ductances defined by slug flow or entrained droplet flow (with poor 
contact between drops and the boiler tube ’vail) . To improve the con- 
ductances by a factor of a hundred or a thousand, rotational flow is 
generated by tube inserts to provide annular flow or high radial 
acceleration droplet flow . 
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BOILING CONDUCTANCE, 



Figure 3-34. Schematic representation of the limits of forced 
flow boiling conductance functions. 
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NOMENCLATURE 



tube cross sectional area 
a numerical constant 
liquid specific heat at constant pressure 
tube inner d meter 

tube equivalent diameter with helical insert 
diameter of wire for wire coil insert 


a function (see Figure 3-17) 
a function (see Figure 3-19) 
a parameter (see Figure 3-26) 



h ... = conductance for an unbroken annular film 
an. film 




boiling heat transfer conductance 


h 


h 


helical flow neat transfer conductance 
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h v (X = l) = 


vap 

k < ■ 

k t - 

k 

v 

L 

I 

Nu = 

N "h * 

Nu < * 
n = 


Pe = 

Pr = 

Pr = 
v 

p 

Re = 


R 'i ■ 


linear flow heat transfer conductance 
heat transfer conductance for a smooth tube 
conductance for dry areas 

conductance for flow at 100 percent vapor quality 
conductance due to vaporizing droplets 
liquid thermal conductivity 

apparent thermal conductivity of turbulent liquid film 

vapor thermal conductivity 

latent heat of vaporization 

tube length 

Nusselt modulus 

Nusselt modulus in helical flow 

Nusselt modulus for linear flow 

number of droplets in tube length of one diameter 
entrained in vapor at reference quality 

Peclet modulus 

Prandtl modulus 

Prandtl modulus of vapor 

absolute pressure, pitch of helical insert 

Reynolds modulus 

Reynolds modulus based on axial flow in the tube 
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Reynolds modulus evaluated with vapor properties 
radius of gas-liquid interface in annular flow 


tube inside radius 


fluid mixed mean absolute temperature 


tube wall absolute temperature 

W | 

superficial liquid velocity, ~ — r 

^ l 

1 

normalized superficial liquid velocity, u^ |gD(p^ - p y )J 2 


superficial vapor velocity. 


normalized superficial vapor velocity, u^|gD(p^ - p y )j 2 
rotating component of velocity at wall (Figure 3-37) 


mass flow rate 


liquid mass flow rate 
vapor mass flow rate 


vapor quality 


vapor quality at tube exit 


a reference vapor quality 


0-571 0.143 


v ^ i i 

Martinelli parameter. Figure 3-5, X = ( — ) \ / ± _ 

tt 'p ' fJ ' 1 x 
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GREEK SYMBOLS 


a 

4p 

4p 


* 



void fraction, contact angle of liquid on wall 


pressure drop 

dimensionless two-phase pressure drop. 


- - p g 

' dz 1 v B 


g°(P £ - P v ) 


- ( 


dP 


dz 


dimensionless single phase liquid pressure drop, gD ^ - p ) 


A P 


s 



* 


A P 


v 


pressure drop in a smooth round tube 
two phase flow pressure drop 


dimensionless single phase vapor pressure drop. 


_ap 


<ZD(P t - p v > 


4p(x = 1) = single phase now pressure drop for 100 percent vapor quality 

4T = wall to saturation temperature difference 
sat 

6 = annular liquid film thickness 

€ = root mean square surface roughness 

£ = Weisbach friction factor 

£ = adiabatic friction factor 

ad 
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«h 


n 


i 




v 


7T 




<t> n 



friction factor for helical flow 

friction factor for a smooth pipe 

absolute liquid viscosity 
absolute vapor viscosity 
3.14159. . . 
liquid density 
vapor density 

a function of tube length to diameter ratio (Figure 3-21) 

a function of wall to bulk absolute temperature ratio 
(Figure 3-22) dP 

( d z) 

Martinet li two-phase pressure drop multiplier, ^p — 

<aF>! 


DERIVATIVES 


M- 

dz 


axial pressure gradient for two phase flow 


d p 

( q 2 ) = axial pressure gradient for liquid flowing alone in tube 
f 

dp 

( "q z ) - axial pressure gradient for vapor flowing alone in tube 
v 
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CHAPTER 4: FLOW REGIME TRANSITIONS 


In Chapter 3, information was provided for defining quantitatively the 
heat transfer and frictional pressure drop of two phase flow in boilers, 
provided that the phase distributions were known. These experimental, 
semi-empirical and analytical data were usually presented in terms of 
single phase functions (for one or both phases). The functions contain 
single phase relationships modified by terms dependent upon vapor 
quality, geometry, flow rate, and other parameters. These results 
demonstrated that the apparently complicated two-phase processes can 
be related to familiar single phase processes. 

In the present chapter, information is provided that will assist the designer 
in determining what phase distributions exist along the length of the 
boiler tube. The phase distribution transitions of principal interest to 
the boiler designer are listed below under five categories, namely, 

(I) fluids which w et the tube wall, (2) fluids which do not wet the tube 
wall, (3) boiler tube inserts, (4) boiler tube thermal conductivity and 
wall thickness, and (5) boiler tube wall surface chemistry and physics. 

The contents of these categories (as they appear in the chapter) are 
outlined below. 

A. Fluids Which Wet the Tube Walls 

1. Superheating of the liquid without phase change 

2. Onset of subcooled boiling 
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3. 


Transition from bubbly flow to wetted-wall slug flow 

4. Transition from wetted-wall slug flow to annular flow 

5. Onset of entrainment from waves on the annular film 

6 . Onset of nucleate boiling in an annular film 

7. "Burnout, " "Boiling Crisis, ■ and "Critical Heat Flux" 

8. Termination of the annular film 

B. Fluids Which Do Not Wet the Tube Walls 

1. Onset of nonwetting slug flow 

2. Onset of transition from slug flow to helical rivulet flow 
C- Boiler Tube Inserts 

1 . Some aspects of two-phase flow in insert transitions 

2. Transition from a helical plug to a helical wire coil 

3. Interrupted inserts 

D. Boiler Tube Thermal Conductivity and Wall Thickness 

E. Boiler Tube Surface Chemistry and Physics 

1. Nucleation sites 

2. Droplet vaporization 


126 



PHASE DISTRIBUTIONS TRANSITIONS 


A . Fluids Which Wet the Tube Walls 

1* Superheating of the liquid without phase change 

The phenomenon of nucleate boiling in conventional pool boil- 
ing systems is attributed to minute cavities with residual gas 
trapped inside. These cavities form the visible nucleation 
sites which are the source of vapor bubbles. Water usually 
boils iii metal containers with only a degree or two of super- 
heat at the wall. However, it has the ability to superheat in 
the liquid phase to almost 150° F above the saturation tempera- 
ture if carefully degassed and heated in a system which has 
been evacuated to remove the trapped gas from the cavities 
in the heated surface. In contrast to metal, glass usually 
has few cavities in the surface which can act as nucleation 
sites. Boiling chips, which contain many microscopic cavities, 
are frequently added to water boiling systems made of glass- 
ware to reduce liquid superheating with its accompanying 
violent change of phase. 

Many boiling systems must be carefully cleaned and evacu- 
ated before they are charged with the working fluid. This 
procedure is common in alkali liquid metal systems, for 
example, and also in systems which use oxidation-sensitive 
organic fluids. The removal of inert or noncondensible gases 
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from the system by evacuation may destroy the effectiveness 
of nucleation sites. 

It is possible to estimate the amount of superheat required 
to bring on the formation of vapor cavities in the liquid in 
the complete absence of nucleation sites. An analysis of a 
model based on the thermodynamic properties of the fluid has 
been made, and predictions presented for water, potassium, 
and sodium in Reference 31. The analysis examines the con- 
ditions for existence of a vapor bubble in a superheated liquid, 
with temperature equilibrium between the vapor and the liquid 
in contact with it. The pressure within the bubble is de- 
termined by the liquid surface tension and the container 
pressure. Similarity between the well-established behavior 
of water and the behavior of the liquid metals is assumed in 
order to predict the behavior of potassium and sodium. These 
results, which represent maximum theoretical values, are 
reproduced in Figure 4-1. 

Many experiments b-.ve been performed to determine the 
amount of superheat required to bring on boiling in forced 
convection systems. These have usually been performed with 
liquid metals. Measured superheats are generally lower than 
those predicted by the analysis (probably at least partly because 
system purity is not perfect), but superheats of hundreds of 
degrees are still common. Experiments performed with 
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potassium (Reference 1) indicate that the amount of superheat 
required to initiate phase change decreases with increasing 
saturation temperature, which agrees with the analytical 
result. 

Apparently a high flow turbulence intensity also inhibits super- 
heating, because these same experiments indicate decreased 
superheat with increasing flow Reynolds number. These same 
trends have been observed in systems containing sodium 
(Reference 2). 

2 . The onset of subcooled boiling 

Since the wall in the preheater is hotter than the bulk mixed 
mean temperature of the fluid, boiling begins at the heated 
surface before the bulk fluid reaches the saturation tempera- 
ture. For ordinary fluids, the mixed mean fluid temperature 
can be significantly subcooled. Initially the vapor is confined 
to a bubble layer adhering to the surface, because bubbles 
separated from c’le wall immediately collapse as a result of 
condensation in subcooled liquid. As the bulk temperature 
increases along th? tube, there comes a point where bubbles 
separated from the wall do not collapse. Downstream of this 
point the void fraction rises rapidly. A diagram illustrating 
this feature is shown in Figure 4-2. In Reference 3, an 
analysis is presented for the liquid bulk subcooling at the onset 
of net vapor generation. The analysis, which is derived for 
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VOID FRACTION 



Figure 4-2. Void fraction variation in subcooled boiling. 
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vertical upward flow, is based upon a balance between the 
forces tending to hold the bubble to the wall and the forces 
tending to remove the bubble from the wall. A contact angle 
function, which describes the surface tension restraining force, 
is evaluated experimentally. Major idealizations include: 
the friction factor of the flow on the bubble layer can be ob- 
tained from classical single phase flow experiments for 
turbulent flow with rough walls, using the bubble radius as 
a roughness element size; Martinelli's analogy solution for 
single phase velocity and temperature profiles holds for the 
bulk liquid flow; and the liquid in contact with the bubbles is 
at the local saturation temperature. 

Comparisons of the results of the analysis with experimental 
data for water and Freon 22 are shown in Figures 4-3 and 
4-4. It may be seen that the comparison is favorable. 

There are several qualifications on the use of this analysis. 

The foremost is that the solution is not appropriate for fluids 
which have very high levels of dissolved gas in them, nor 
for highly degassed systems, since both of those will exhibit 
abnormal nucleation behavior. Another qualification is that 
the boiling surface must have a wide size range of active 
nucleation sites available. 

In Reference 4, relationships are derived for the vapor void 
fraction in subcooled boiling of ordinary fluids, utilizing an 
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Figure 4-3. Comparison of subcooled boiling water data with analytical 
predictions (adapted from reference 3). 
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Figure *-•*. Comparison of subcooled bailing Freon 22 data with 
analytical predictions (adapted from Reference 3). 
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analytical fit to experimental data for the mixed mean tempera- 
ture distribution. The temperature distribution function is 
shown in Figure 4-5. With this temperature distribution and 
the known heat flux, the void fraction is obtained from an 
energy equation and the initial point of net vapor formation. 

The relationship for void fraction is compared to experimental 
data, both for negligible subcooling conditions, where there 
are many sets of data, and for cases with significant subcooling 
as well. A typicai comparison is shown in Figure 4-6. It 
may be seen that the analysis agrees quite well with the data. 

The quantity C , given in the figure, is an experimentally 
o 

determined phase distribution parameter. 

3. Transition from bubbly flow to wetted-wall slug flow 

Bubbly flow, a two-phase mixture of many discrete small 
bubbles in a primarily liquid flow, is easily produced in 
adiabatic systems. As the void fraction increases, this phase 
distribution gives way to one in which many of the bubbles 
coalesce to form vapor slugs separated by bubbly-liquid bridges. 

A number of investigators have produced maps of flow 
regimes in adiabatic systems which include this transition. 
However, in actual boiling systems utilizing ordinary fluids, 
the bubbly phase distribution apparently exists only in the sub- 
cooled soiling regime (Reference 6). At inception of vapor 
quality, .he bubbly flow is interspersed with vapor slugs, so 
that in the phase distribution descriptions used in this report, 
it would be referred to as a wetted-wall slug flow. 
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Figure 4-6. Comparison of experimental and predicted void fraction for 

subcooled water boiling, C 1.13. 
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In Reference 7, the authors conclude, on the basis of a 
literature survey and their own experiments, that for adiabatic 
flow there is no simple criterion which can be used to charac- 
terize this transition, but that wetted-wall slug flow always 
occurs at void fractions above 35 percent. In many boiling 
systems utilizing ordinary fluids, a void fraction this high 
occurs well before saturated boiling starts (see, for example. 

Reference 4)- It appears that the vapor generator designer 
may assume that this transition occurs near onset of satu- 
rated boiling without serious error in either pressure drop 
or heat transfer conductance calculations . 

4. Transition from wetted-wall slug flow to annular flow 

This transition, which occurs at low' vapor quality, has been 

studied anal lically and experimentally by a number of labora- 

8 9 10 11 

tories, usually utilizing adiabatic systems. ’ ’ ’ The 

transition is not a particularly sharp one, and several sources 
describe other phase distributions in between. 

A flow regime between slug and annular flow frequently de- 
scribed is "churn" or "semi-annular" flow, which is a very 
unsteady oscillation between annular flow and slug flow in 
adiabatic systems . However, since at the higher quality end 
of wetted slug flow the slugs become longer and longer so that 
the major portion of the flow is indistinguishable from annular 
flow from a heat transfer point of view, relationships which 
describe this transition only approximately are probably 
adequate. 
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The analysis of the transition from slug to annular flow has 
been undertaken by several authors, using the approach of 
finding the intersection of analytical expressions for void 
traction of both the slug flow and annular flow regimes. Because 
of the difficulty in representing either of these regions ana- 
lytically, and because the void fraction at the onset of annular 
flow is very little effected by fluid properties and flow rates 
(a is close to 0.80 for onset of annular flow for most ordinary 
fluids), these attempts at a theoretical approach have not 
been particularly successful. 

A large body of adiabatic data was accumulated and presented 
in Reference 8. These data were obtained from systems with 
a variety of fluids and tube diameters, using as a criterion for 
the transition the presence or absence of liquid bridges in the 
flow. 

On the basis of these data, two correlations for the slug- 
annular transition are presented in Reference 8. These 
correlations are shown in Figure 4-7. The sharp break in 
slope is attributed to the onset of spray entrainment in the 
vapor flow. Also reported in Reference 8 are data on void 
fraction and pressure gradient at transition for some specific 
systems (Figure 4-8). The pressure gradient goes through 
a rather indistinct minimum in the region of this transition. 

In Reference 11, data for upward adiabatic flow of steam- 
water mixtures are reported. These data have been used 
to construct a graph of the slug-annular transition specifically 
for water-steam mixtures as a function of saturation pressure. 
This graph is shown in Figure 4-9. 
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Figure -4-7. Correlation for transition from slug to annular flow. 
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Figure 1-8. Void fraction and pressure drop at transition from 
slug to annular flow. 
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Figure 4-9. The slug-annular flow transition for steam water 
mixtures in upward flow. 
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5. Onset of entrainment from waves on the annular film 

One of the important mechanisms causing the entrainment 

of liquid into the vapor phase from the liquid layer in annular 

flow boiling is vapor shear. It appears that most of the shear 

entrainment is from the peaks of large "breaking" waves, 

which have a height above the average liquid layer thickness 

12 

of the same order as the layer thickness itself. There is 
no agreement as to the conditions upon which the onset of 
roll waves occurs. However, there has been a correlation 
published for the onset of spray entrainment from these waves 
for adiabatic air-water systems. 

A simple analysis, presented in Reference 3, allows the 
prediction of shear entrainment from other fluid systems, 
and agrees quite w ell with the correlation referred to above 
for air-water. It is based on the following model. 

Presume that the waves already exist on the liquid surface, 
and that the flow is upward. The removal of a droplet is 
supposed to be due to the drag of the vapor upon the tip of 
an irregularity in the wave and the restraining forces con- 
sidered are surface tension and gravity. A schematic drawing 
of this model is shown in Figure 4- 10a. The force system 
acting upon the wave tip is indicated in Figure 4-10b. Clearly 
the restraining force on the wave tip decreases and the drag 
force increases rapidly after a "neck" appears, and a droplet 
is almost certain to be pulled off. 
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Figure 4-10. Shear entrainment model for liquid and vapor 
in upward flow. 
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Exaluation of a force balance at the point in time that the neck 
is beginning u form, utilizing the simplifications that the 
liquid velocity is negligible compared to the vapor velocity, 
that the resulting droplet diameter is equal to the mean liquid 
layer thickness, and that the drag coefficient for the downstream 
end of a rounded end cylindrical object is a constant, leads 
to a closed form relationship for the critical vapor velocity. 

This relationship is given in Figure 4-11. 

By combining this shear entrainment criterion with the annular 
flow liquid film thickness predictions from Chapter 3, a map 
of conditions under which shear entrainment would be expected 
results. Two such maps are presented in Figures 4-11 and 
4-12. 

As an example of the interpretation of these figures, consider 
the graph for Freon 12 flowing in a 11.5 mm tube, Figure 
4-11. By hypothesis, the film thickness must be at least 
equal to the droplet diameter (or wave height) for entrainment 
to occur. Therefore, a given size droplet will only be en- 
trained into the vapor at qualities for which the film thickness 
is greater than that drop size. However, the vapor velocity 
must be high enough to entra droplets of that size so the 
entrainment region lies above the line representing the en- 
trainment criterion. Conditions for entrainment of 1.3 mm 
Freon droplets exist above the line of the shear entrainment 
criterion, and below the 1.3 mm film thickness line, in the 
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Figure* 4-11. Boundaries for onset of entrainment due to shear. 
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Figure -*-12. Boundaries for onset of entrainment due to shear. 
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cross-hatched area. Figure 4-12, which is constructed for 
wetting flow mercury, indicates that shear entrainment is 
extremely unlikely under conditions of practical importance. 

6. Onset of nucleate boiling in an annular film 

This transition is important for several reasons. One of 
these reasons is that nucleate boiling alters the heat transfer 
conductance associated with annular flow, as shown in Chapter 
3. A second reason is that nucleate boiling ejects liquid from 
the film and, therefore, increases the entrained liquid. The 
loss in liquid from the film may lower the quality at which 
the annular film ends. 

From the models of heat transfer in annular flow given in 
Chapter 3, one would expect that a thick annular film would 
exhibit nucleate boiling at a lower heat flux than would a thin 
one: i.e. , nucleate boiling would appear at higher flux as 
the quality increases. One would also expect that a low- 
turbulence film would exhibit nucleate boiling at a lower heat 
flux than would a very turbulent one; i.e. , at fixed quality, 
a higher flow would require higher flux to bring about nu- 
cleation. These trends both appear to be corroborated by 
the data ir Figure 4-13, which are for water-steam at 7 atm 
taken from Reference 14. The test section for the experi- 
ments aa ~ heated rod surrounded by a glass tube to form an 
annular flow passage. 
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If information is available concerning the amount of wall 
superheat required to bring about nucleate boiling, then by 
use of the methods of Chapter 3 the heat flux at which nucleate 
boiling will occur can be predicted and, in the complete absence 
of nucleation sites, the liquid superheat predictions of Section 
A. 2 of this chapter are applicable. 

As an example, consider the conditions under which nucleation 
will start in an annular flow r of a liquid metal. In the absence 
of nucleate boiling (which acts as a volume heat sink within 
the liquid film), all of the heat added to the annular film at 
the wall must be conducted or convected across the liquid 
layer to vaporize fluid at the inner surface, which is taken 
to be at the saturation temperature. Since the liquid metals 
have very high thermal conductivity, convection may be neglected 
and the temperature drop across the liquid layer computed on 
the basis of conduction alone. Nucleation will occur when the 
heat flux is adequate to maintain the liquid in contact with the 
wall at the temperature for nucleation given in Section A. 2. 

With this simple model, a graph for the minimum heat flux 
required to bring on nucleate boiling has been computed for 
potassium and is shown in Figure 4-14. A similar curve has 
been constructed for Freon 12, using an arbitrary superheat 
temperature. The curve for Freon is, of course, hypothetical, 
since with its low thermal conductivity, turbulent convection 
in the liquid layer certainly predominates (see Figure 4-15). 
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Figure 4-14. Heat flux required to initiate nucleate boiling in an 
annular flow without nucleation sites. 
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Figure 4-15. Heat flux required to initiate nucleate boiling in an 
annular flow without nucleation sites. 
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7. "Burnout, * "Boiling Crisis, " or "Critical Heat Flux" 

If, in an experiment, a boiler with fixed geometry and fixed 
mass flow rate is subjected to an increasing heat flux, an 
operating point will be reached at which, for a very small heat 
flux increase, a large increase inwall temperature will occur. 
Usually, the wall temperature will not increase suddenly over 
the entire boiler. If a stable operating point is obtained at this 
new condition, it will usually be found that at some axial location 
near the outlet end the wall temperature will change, over a 
relatively short distance, from a low value of a few tens of 
degrees to a high value of a few hundreds of degrees. In a 
system with a uniform heat flux, the boiler tube may fail at 
this point. There have been many specific experiments 
investigating this phenomenon, which is commonly referred 
to as "burnout, " "critical heat flux," or, less commonly, 
"boiling crisis. " 

Unfortunately, the notion of a "critical heat flux" can lead 
to some erroneous concepts concerning the phenomena being 
observed. Several phenomena are actually represented. One 
of these, film boiling, is described in most heat transfer texts. 
Factors controlling its onset are described in Section E of this 
chapter. However, probably the most important of the 
phenomena associated with "burnout" are dryout of the annular 
liquid layer on the wall or breakup of the liquid film into rivu- 
lets. Of course, a forced convection boiler designed to produce 
dry vapor at the exit must operate through this region without 
failure, so the term "burnout" is inappropriate. 
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Observations made with Freon 11 in a glass test section 
(Reference 15) support the interpretation of "burnout” as 
a breakdown of the annular film. These observations, taken 
under uniform heat flux conditions, demonstrate a very large 
and discontinuous decrease in heat transfer conductance at an 
axial location which coincides with the dryout of the annular 
film. 

The rate of change of the liquid flow in the film is determined 
by a mass balance between the liquid evaporated, the liquid 
deposited from the vapor, and the liquid entrained from the 
film. The liquid entrained from the film is usually supposed 
to be removed from the film by stripping of liquid from wave 
peaks, and ejection of liquid by breaking of bubbles generated 
by nucleation in the film. The mechanisms controlling the 
flow rate in the annular film are, therefore, largely hydrodynamic. 
One might expect the liquid film flow rate at any point to be de- 
termined by the processes occurring upstream; the rates of 
deposition . id entrainment, and the rate of change of vapor 
velocity all contribute to the liquid film flow. Many different 
combinations of these processes would give virtually the same 
liquid flow rate at a particular point. The annular film is very 
thin near its end so that bubble nucleation and wave motion are 
probably both suppressed, and the principal mass fluxes in and 
out of the annular film near its termination are probably due to 
evaporation and deposition. 
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When film dryout occurs at a quality less than unity, 
the remainder of the liquid mu. - -t certainly be entrained in 
the vapor flow. Due to diffusion, the fine droplets in the vapor 
eventually settle out on the walls and evaporate. The longer 
the boiler, the smaller the entrained flow at the end will be. 

It is clear that for a fixed flow rate, a low heat flux boiler is 
a long boiler. If entrainment rates are not too high, a long 
boiler has proportionately less entrained liquid at a given quality, 
and one would expect that the end of the annular film would 
occur at higher quality for a low flux (long) boiler than for a 
high flux (short) boiler. Regardless of tube diameter, liquid 
flow rate, or other conditions, the end of the annular film will 
occur at X = 1 for vanishingly small heat flux. 

The plot of Figure 4-16 (Reference 13), showing Freon 12 
burnout data, indicates a linear relationship between (q/A) 

B 

and X over a wide range of test conditions. However, all 
B 

of the data lines must pass through the point (q/A) = 0, 

B 

X = 1. Since this is the case, only one of the lines in 
B 

Figure 4-16 can be straight at substantially lower flux than 
the data cover. All of the remainder must curve at low heat 
flux. The data of Figure 4-16, plotted on semi-logarithmic 
coordinates, are shown in Figure 4-17. It is clearer in this 
plot that the data do tend toward zero heat flux at X = 1. The 
solid lines indicate the r.ctual range of the data. The dotted 
lines indicate where one might expect data taken over a wider 
range of he^t fluxes to lie. These same data have been replotted 
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Figure 4-16. Freon 12 linear flow "burnout" data (from reference 13). 
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as functions of length in Figure 4-18. It may be seen that the 
boiler length required to maintain the annular film at high quality 
increases rapidly for high mass flow rates. 

As stated earlier, when the wall is no longer wet, all of the 
remaining liquid is entrained in the vapor. An entrained liquid 
plot (Figure 4-19) has been constructed from the data of Figure 
4-18. It can be seen that for high flows and fluxes, as much as 
90 percent of the liquid can be entrained 

8. Breakdown of the annular film 

Stability Analyses 

A number of authors have studied the breakdown of liquid 
films flowing over heat transfer surfaces. In References 
16 and 17, the authors studied the flow of water downward 
over vertical heated surfaces. They concluded that the de- 
velopment of longitudinal striations (described as "roping") 
was due to the variation in surface tension with film thickness. 

In Reference 18, the author treated analytically the stability 
of downward flowing films subjected to heat transfer on the 
basis that the destabilizing mechanism is the nonuniform 
pressure normal to the surface due to variations in evaporation. 

The analysis led to definition of a critical heat flux for destabi- 
lization, and a prediction of the rate at which dry spots will 
form after the breakdown process begins. 

Although these works are not strictly applicable to forced convection 
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Figure 4-19. Entrained liquid at end of annular film - Freon 12. 



















boiling since the interface shear present in forced convection 
boiling does not appear, extension of these approaches to an 
appropriate case is probably feasible. 

Other authors have studied the stability of dry patches on 

19,20,21,22 

walls in annular flow. An analytical treatment of 

film stability on unheated walls, given in Reference 21, was 
compared to adiabatic experimental data in Reference 20. The 
experiments were performed by blowing an air jet on the 
annular film to expose the wall and then deterr.- ning the con- 
ditions under which the wall would rewet when the air jet was 
removed. The rewetting depsnds strongly upon contact angle 
uclween the dry wall and the liquid- The comparison indicator 
that aerodynamic forces on the edge of the liquid are import 
In Reference 19, which treats the dry spot under conditi -f 
heat transfer, it is shown that for liquids ith high wettability, 
the thermal effects become dominant in determining the sta- 
bility of dry patches and that surface forces dominate in cases 
of poor wettability. In Reference 5, it is shown that thermal 
conduction in the wall can be neglected only under conditions 
in which wall temperature gradients near the boundrry of the 
dry spot are small. 

Another approach to the onset of rivulet flow is presented 
in Reference 5 . For this analysis, it is postulated that the 
vaporizing liquid is in direct contact with the smooth tube 
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surface and the boiling consists of surface vaporization or 
some volume boiling but no film boiling If the contact 
angle between the liquid and the wall is zero degrees and 
there is no entrained liquid then, hypothetically, there will 
be a complete annular liquid film at the wall from qualities 
ranging from zero to 100 percent. If the contact angle (which 
is defined by the liquid-vapor, solid-vapor, solid-liquid surface 
forces) is greater than zero, the width of a liquid rivulet for 
that particular profile area defined by the contact angle and 
vapor quality must be compared to the inner circumference 
of the boiler tube If the rivulet width is greater than the tube 
circumference, then the liquid film would not rupture. If the 
rivulet width is less than the tube circumference, then the 
liquid film would be separated. 

Rough tube surfaces can affect the liquid film in several ways. 
Surface asperities which are high enough to penetrate the film 
free surface reduce the liquid film contact area, as do non- 
hon. .geneous structures such as oxide deposits, which change 
the local contact angle. 

Finally, film boiling can occur at surface pits or local areas 
of nonwetting, reducing the wall superheat temperature differ- 
ence at which surface vaporization ends ami transition boiling 
begins . 

In Reference 5 , an analysis is presented for the conditions 
•■.nder which rivulets will form on a smooth tube surface- 
The critical contact angle predicted in this analysis is only 
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a function of tube radius and thickness of a complete annular 
liquid film (see Figure 4-20). From the relation between 
vapor quality and the annular film thickness, it can be shown 
that the critical contact angle decreases as vapor quality 
increases. Thus, one would expect the annular film to 
rupture at vapor qualities less than 100 percent. A specific 
example showing this effect is given in Figure 4-21. 

Film Flow Rate Analyses 

In Reference 23, the authors present an analysis of the film 
flow rate in the presence of vaporization, entrainment and 
deposition in an annular-dispersed flow. The analysis, 
which treats the heat transfer term as a per*urbation on the 
adiabatic case, depends upon experimentally-determined 
adiabatic data for its evaluation. The deposition rate is 
treated analytically, utilizing the turbulent diffusivity 
from single phase turbulent flow relationships to establish 
the deposition flux. Because of the nature of the data with 
which their comparison is made, the authors use experi- 
mentally-determined film end points to predict "burnout 
heat flux" and then compare that value to the experimental 
flux. The comparison with selected boiling water data for 
a variety of heat flux distributions shows that trends are 
correct, even though the analytical approach is simple. 
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Figure 4-21 . Critical contact angle versus vapor quality for a 
specific potassium system. 
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B. Fluids Which Do Not Wet the Tube Walls 


1 . Onset of nonwetting slug flow 

Fool boiling mercury in an unwetted glass vessel at re- 
duced pressure exhibits a behavior which is in marked 
contrast to boiling water, a wetting fluid, under the same 
conditions . When the mercury reaches the point of incipient 
boiling, a vapor bubble appears against the glass, then rapidly 
expands into a thin sheet of vapor against the wall. This sheet 
of vapor then slides upward along the v;all to the surface. It 
is clear that, under these circumstances, the path of least 

resistance for the vapor is not through the liquid as a train 

24 

of bubbles, but along the wall. 

Although some of the properties of mercury are somewhat 
extreme (very high surface tension and density), and the 
conditions under which the observations were made are 
unusual, the phenomenon observed must be largely dependent 
upon the wetting properties of the fluid-wall combination. A 
somewhat similar "film boiling" may be visualized for the 
onset of boiling in any nonwetting flow. When vapor is form- 
ed, the liquid may be pushed away from a significant portion 
of the heat transfer surface, so that the flow tends toward 
agglomerated masses of liquid separated by volumes of 
vapor. This is the distribution which characterizes non- 
wetting slug flow. A designer working with nonwetting 
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fluid-wall combinations can probably conservatively assume 
that the initiation of slug flow coincides with the onset of 
phase change and experiments in small diameter tubes with 
mercury flow confirm this (Reference 15). 

2. Transition from slug flow to helical rivulet or droplet flow 
Heat transfer conductances for mercury boiling in a tube 
containing a helical divider are shown in Figure 4-22, taken 
from Reference 15. It may be seen that at low quality the 
heat transfer conductances have the same general shape as 
those for slug flow (Figures 3-7 and 3-8). As the quality and 
vapor velocity increase, the conductances undergo a change, 
increasing by an order of magnitude or more. This change 
is associated with the onset of helical rivulet or droplet flow. 

At low vapor velocity, the liquid velocity is very low, and 
little radial acceleration is provided by the helical divider. 
As the liouid velocity increases, it spreads out into a film 
boiling rivulet which is held against the wall by the radial 
accelerations. The heat transfer conductances after the 
transition may be predicted from the model of Chapter 3. 

The transition is, of course, gradual, since droplets or 
globules of different sizes have different amounts of slip 
with regard to the vapor velocity. Experiments in visuali- 
zation of boiling mercury (Reference 25) have shown 
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Figure 4-22. Helical flow boiling mercury conductances versus quality. 
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that range of globule-droplet sizes in upward flow of 
mercury is quite large. The criterion for the beginning 
of this transition is that the radial acceleration for the 
smaller droplets, including slip, must be many times the 
acceleration of gravity. 

C. Boiler Tube Inserts 

1. Some aspects of two phase flow in insert transitions 

Axial changes in boiler insert type which have been con- 
sidered for use by designers include: fluted plug to helical 
vane, fluted plug to helical wire, and helical vane to helical 
wire. The vane and wire sections may or may not have a 
centerbody associated with them, and the insert pitch may 
or may not be the same on both sides of the insert change. 

Consider the transition from a fluted plug to a helical vane 
system. If tne pitches of both are the same, then the major 
difference between the two inserts is in the diameter of the 
centerbody, and the flow will be similar to a flow into a sudden 
enlargement. The vapor phase will decelerate rapidly, and 
there will be a recirculating stall at the base of the plug if the 
area change is large. The liquid layer on the wall will also 
tend to lose velocity and thicken, although at a lesser rate 
than the vapor, and if the quality change with length is high 


169 



enough, no significant thickening of the liquid layer will occur. 
No observable effect of this transition upon the heat transfer 
performance would be expected. 

If the fluted plug and helical vane do not have the same pitch, 
the liquid layer will be subjected to an acceleration which will 
thicken it on one side of the channel formed by the next insert 
and thin it on the other side. The liquid layer flow is similar 
to the open channel flow around a sharp corner, except that 
the gravitational field is replaced by a centrifugal acceleration 
field which may be many times larger than gravity. From 
experience with open channel flows, one would expect that the 
major effects of the sharp corner turn would be obliterated 
within a couple of channel widths. 

One would expect a liquid layer on the tube wall to have a 
deleterious effect upon the flow rotating effectiveness of 
only the helical wire. If the liquid forms a layer which has 
appreciable thickness when compared to the diameter of the 
wire then, as was shown in Chapter 3, the force tending 
to rotate the fluid will be decreased, and if the liquid layer 
covers the wires, no force will exist at all. However, the heat 
transfer downstream from a transition to an ineffective wire 
coil will not necessarily be poorer than that in a helical vane. 
As long as the wall is completely covered, liquid-like con- 
ductances will prevail. Since centrifugal accelerations are 
decreasing, more liquid may be entrained in the vapor, and 
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surface waves will also be deeper. Both of these effects 
thin the liquid layer and increase turbulent diffusion. If 
nucleation does not occur, then the vaporization is at the 
liquid-vapor interface, so that a thin turbulent layer, which 
has less heat transfer resistance than a thick one, will 
exhibit a higher heat transfer e inductance. Decreased 
radial acceleration also decreases the pressure rise across 
the film, which can inhibit nucleation. 

As long as the liquid film on the wall remains intact, none 
of the insert transitions would be expected to cause large 
changes in the heat transfer or pressure drop characteristics 
of the boiler. 

2. Transition from a helical plug to a helical wire coil 

’-1 Chapter 3, relationships were presented for the fully- 
established flow of a single phase fluid in a helical wire coil. 

It is, however, under relatively unusual conditions that fully- 
established flow will be obtained because of the very long 
entrance lengths required. In Reference 5, an analysis was 
presented for transition to fully-established flow in a wire 
coil, starting from an arbitrary initial condition. The model 
is one of unsteady solid body rotation, with forces exerted 
on the (solid) fluid cylinder by the wire drag and wall friction. 

The analysis leads to relationships for the asymptotic approach 
to fully-estabUshed flow. In order to illustrate the nature of 
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the results, a specific case is presented in Reference 5. The 
axial distribution of rotational velocity for this case, which is 
transition from a helical vane to a helical wire coil, both of 
which have the same pitch, is shown in Figure 4-23. The 
characteristic length for this transition is shown in Figure 
4-24. As may be seen in the figure, many diameters of tube 
are necessary for flow accommodation. 

3. Interrupted inserts 

Many experiments have t.hrwn that wetting liquids tend to 
collect on the centerbodies of helical vane and multi-groove 
inserts, because of the secondary flows . An analytical 
investigation of this effect which corroborates experimental 
results, is presented in Reference 5. Because the center- 
body is usuallv not a heat transfer surface, the accumulated 
liquid does not readily evaporate, and can run along the 
centerbody far into the superheat region. For this reason, 
many designers who have utilized vane systems with center- 
bodies have interrupted these inserts at intervals in order to 
return the liquid to the wall. When the liquid comes to the 
end of the centerbody it is stripped off into the rotating vapor. 

In order to predict the necessary spacing for these interruptions 
in the insert, an analysis of the droplet motion has been worked 
out in Reference 5. 
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Figure 4-23. Change in rotational velocity of a flow passirg from a helical vane 
insert to a helical wire insert, both of the same pitch. 
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Figure 4-24. Normalized characteristic length. 
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Several phenomena are taking place simultaneously which affect 
the droplet motion. The droplets, im .. .'.ly at res';, slip with 
respect to the vapor in both the axi;> mil rotational motions, 
and are accelerated in a three-dimensional sp<r-.l path due to 
drag forces . At the same time, the vapor rc >n, no longer 
maintained by the helix, begins to decay due to boiler wall 
friction, and the vapor axial velocity continues to increase due 
to vapor addition from vaporization. 

In addition, the droplets may be undergoing fragmentation due 
to internal flow instabilities, vapor shear gradients, or other 
fcioes. In order to analyze this complex two-phase flow, a 
number of restrictions have been made, which are: 

(a) The axial vapor velocity field is uniform across the 
tube, and its change is negligible witl n the axial 
length . 

(b) The vapor angular velocity change is negligible within 
the axial length required for the droplets to impact 
on the wall. 

(c) The droplet ~ize is invariant during flight. 

(d) The vapor rotational velocity is a linear function of 
radius (solid body rotation). 

(e) The droplet drag coefficient is constant. 
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A sketch of the geometry of the composite problem is shown 
in Figure 4-25. 

Two limiting cases have been analyzed. In one, appropriate 
to large droplets, the idealization is made that the droplet 
velocity is always small compared to the vapor velocity, i.e. , 
the slip is large. In the other case, appropriate to small 
droplets, the idealization is made that the droplet slip is at 
all times small. The range of applicability of each of these 
two solutions has been estimated. 

The results of these two analyses are shown in Figure 4-26. 

It 1 . .ay be seen that the two solutions provide a reasonable 
bound for the axial travel distance of entrained droplets. 

The analytical result for the small droplets is given in the 
Appendix. The result for the large droplets is not in closed 
farm. 

A rather striking photograph of the deposition of droplets 
shed from the interrupted centerbody of a helical insert is 
shown in Figure 4-27. This photograph, which suppotts the 
analysis described above, was taken during the air-water 
experiments reported in Reference 26. 

D. Boiler Tube Thermal Conductivity and Wall Thickness 

Many experiments have shown that the onset of film boiling occurs 
when the wall surface temperature exceeds me liquid saturation 
temperature by a cri teal value. 
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Figure 4-25. Idealized droplet motion at end of centerbody in 
helical flow. 
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Figure 4-26. Axial distance to impact on wall for droplets In helical flow. 



Drop lot deposition downstream from an interrupted 
centerbody (photo courtesy of G. E. Missile and 
Space Division, Cincinnati). 



it was shown in Reference 5 that local heat fluxes and wall tempera- 
tures are much higher than the average values near discontinuities 
in the inner wall conductances. The reason for this is that a dis- 
continuous change in conductance is accompanied by a change in 
inner wall temperature, so that a significant amount of heat flows 
in the wall from the region of low conductance to the region of high 
conductance. Figure 4-28 taken from reference 5 indicates this 
effect for the case of uniform outside wall heat flux. When the 
surface temperature under the edge of the liquid film rises to a 
critical level, local " vapor blanketing 0 occurs (and the end of the 
liquid film moves upstream). 

The closer to the beginning of the boiler the film ends, the higher 
the droplet concentration in the vapor and the deposition rate onto 
the wail must be- The impinging droplets may not wet the wall, 
but they still contribute to the wall cooling, so that the ratio of the 
'wetted wall' heat transfer conductance to the 'dry wall' conductance 
decreases as the end of the film moves toward the boiler inlet. As 
is shown in Figure 4-28, the decrease in this ratio causes a lowering 
of the wall temperature at the edge of the film, so that at some point 
the edge becomes stable and ceases to move. Given a fixed value 
of wall to saturated fluid t<- ..perature difference for which the onset 
of film boiling occurs, the maximur. heat flux for which the end of 
the annular liquid layer is stable can be computed from the data in 
the figure and the definition of the heat transfer conductance. 
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Figure 4-28. Normalized temperature at a conductance discontinuity. 



One conclusion that may be drawn from this line of reasoning is that 
a t'.iick walled or high conductivity boiler decreases the likelihood 
that ti - fa .nular liquid layer will be removed by film boiling. 

*iK..ne-- conclusion is that internal flow regime changes occurring 
over short distances in a boiler tube would appear to extend over 
much greater distances if inside wall conductances are inferred 
from outside wall temperatures. An example of such an interpre- 
tation may possibly have occurred in the boiling Freon data of 
Reference 27. Those data, obtained with a (high conductivity) 
copper tube, show a change from the maximum near the end of the 
annular film to the minimum in the fog flow region which is quite 
gradual, while the data for Freon given in Reference 28, obtained 
with the aid of flow visualization in a (low conductivity) glass test 
section show a discontinuous transition between the two regimes. A 
comparison of these data is shown in Figure 4-29. 

E. Boiler Tube Surface Chemistry and Physics 
1. Nucleation sites 

Several choices which the designer makes strongly effect the 
availability of nucleation sites in the boiler. Not only are the 
wall material and tube wall fabrication techniques important, 
but the methods of cleaning the system and introducing the 
working fluid affect the nucleation performance. Very small 
reentrant cavities filled with gas other than the working fluid 
vapor to form nucleation sites are usually found in profusion 
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NORMALIZED HEAT TRANSFER CONDUCTANCE 



Figure 4-29. Linear forced convection boiling data for Freon 1 and 22. 
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on normal wrought materials. However, surface preparation 
can strongly affect them. Figures 4-30 and 4-31 (Reference 29) 
indicate the magnitudes of the effects which can occur with pool 
boiling liquid nitrogen. In Reference 39, it was shown that 
nucleate boiling conductances of Freon 113 on copper could be 
improved as much as fivefold by surface treatment such as 
sandblasting, electroetching, etc. It would be expected that 
some small improvement in liquid metal boiling could also be 
obtained by such treatments. However, for low Prandtl 
number fluids, the controlling thermal resistance is usually 
not in the boundary layer. 

Surface material effects are also observed in forced convection 
boilers. If nucleate boiling does not occur, then heat which is 
added to an annular film at the wall must be conducted across 
the liquid layer so that vaporization occurs at the interface 
between the liquid and the vapor. One would expect that in the 
absence of nucleation sites, the onset of nucleate boiling would 
be dependent upon the maximum superheat temperature which 
the fluid against the wall can sustain w ithout changing phase. 

In the forced convection experiments of Reference 30, Freon 22 
was vaporized in three tubes of identical geometry and under 
very nearly identical conditions, except that two of the test 
sections were metal and the third was glass. The 0 lass test 
section was heated by electrical conduction through a transparent 
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Figure 4-30. Effect of materials on nucleate boiling curve for 
nitrogen (Reference 29). 
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"igure 4-31. Effect of number of drilled cavities on nucleate-boiling 
cuive for nitrogen (Reference 29). 
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metal film on the inside and was used for flow visualization. 

Some of the data from these experiments are reproduced in 
Figures 4-32 and 4-33. In these data the glass exhibited a 
much lower boiling heat transfer conductance than the two 
metal tubes. Since the glass test section does not have as 
many nucleation sites as the metal tube, more superheat of 
the liquid at the wall is required to bring on vapor generation 
in the liquid bulk. 

The designer can, by choice of surface material, fabrication 
method, and system filling technique, either enhance or 
suppress nucleate boiling in forced convection. 

2. Droplet vaporisation 

The chemical and physical nature of the boiler tube surface can 
control the thermal behavior of vaporizing drops. Therefore, a 
brief review of the mechanisms in operation and the corresponding 
transition criteria are presented.* 

The lifetime of a droplet on a heated surface, i.e. , the length 
of time from its placement on the surface until it has completely 


Many contributions in the field of droplet vaporization can be found in the 
literature (see References 32 - 38, for example). 
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Figure 4-32. 


Heat transfer data fcr refrigerant 22 vapor qualities from 
0.01 to 0.60. 
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4-33. Heat transfer data for refrigerant 22 vapor qualities trom 
0.01 to 0.60. 
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disappeared, may be divided into three regions, depending upon 
temperature difference. At low temperature differences, the 
droplet is in contact with the heated surface, and vaporizes at 
the free surface. As the temperature difference is increased, 
the lifetime, which is proportional to the reciprocal of the heat 
flux, decreases, until it reaches a minimum (end of the surface 
vaporization region) . As the temperature difference is increased 
beyond this point, the transition region begins in which the 
droplet is partially supported on a layer of its own vapor. At 
higher temperature differences, the droplet floats upon a film 
of its own vapoi (like a ground effects machine) with the heat 
transfer taking place by convection and radiation across the 
vapor film. 

Analyses which describe the direct contact surface vaporization 
process are given in Reference 35. The lifetime prediction from 
one of these models describes the vaporization of a droplet 
idealized as flat disk vaporizing from the upper surface only. 
Idealizations required in the analysis are that the upper surface 
is at the saturation temperature, and that the droplet disk 
diameter does not change during vaporization. The disk simply 
grows thinner with time. An interface thermal resistance between 
the droplet and the surface, interpreted as the portion of the disk 
area which is in direct contact with the surface, is included in 
the analysis. 
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If the surface is not microscopically smooth, tbe pits in the 
surface will probably not be completely filled with liquid, but 
will contain a small gas volume. An equation can be written 
relating the pressure and liquid surface tension forces, which 
defines the conditions under which the liquid-vapor interface 
will no longer remain within the pits. Such a relationship is 
presented in Reference 35, based on the approximation that the 
liquid vapor pressure is a linear function of the temperature in 
the range of interest. 

Analyses of the film boiling of droplets above heated surfaces 
under several conditions are also available. In Reference 32, 
the lifetime of film boiling droplets on smooth surfaces is 
presented. The film boiling of a droplet on a rough surface, 
in which some heat is transferred by conduction through 
roughness elements on the surface, as well as convection 
through the vapor film, is presented in Reference 34. 

Under some conditions, the droplet will spall or fission into 
several fragments. The lifetimes of the individual fragments 
differ from the lifetime of the complete droplet. Under some 
conditions, the droplets oscillate in the film boiling region. 

It has been shown analytically in Reference 34 that this motion 
should increase the droplet lifetime. 
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In Figure 4-34 are sho.vn representative analytical curves for 
the models describing direct contact surface vaporization on 
smooth and rough surfaces, film boiling on smooth and rough 
surfaces, spalling and oscillating droplets. 

The presence of dissolved impurities in droplets can change 

the behavior of droplet vaporization in two ways, depending 

upon whether the impurities are surface active or inactive. 

hi the first case, the accumulation of surface active impurities 

(e.g. , potassium in mercury) on the droplet surface reduces 

the vapor pressure or percentage of vaporizing liquid exposed 

to the heat transfer surface relative to a pure droplet, therefore 

35 

changing the heat transfer to the droplet. The presence of 
the surface active impurities effectively increases the lifetime 
of the droplet, hi the second case, the accumulation of dis- 
solved surface inactive impurities (e g., iron in mercury) 
upon the droplet surface as a result of impurity saturation in 
the vaporizing droplet can be visualized to increase the heat 
transfer. The surface inactive impurity does not accumulate 
upon the surface until the droplet becomes saturated. When 
the droplet becomes saturated, the concentration of mercury 
on the droplet surface is suddenly reduced, thereby causing the 
droplet to contact the heat transfer surface. Thus, a change 
in the droplet vaporization mechanism from film boiling to a 


192 




Figure 4-34 . Some of the important boiling heat transfer states for droplet vaporization. 



direct contact heat transfer type results. (This mode of 
droplet behavior has been observed in mercury vaporization 
experiments reported in Reference 35.) 

Experimental data for the vaporization lifetime of a variety 
of different liquids are shown in Figures 4-35, 4-36, and 4-37. 

In some cases, the theoretical prediction for film boiling on 
a smooth surface is also included. A series of droplet vaporization 
experiments is described in Reference 36, to determine the 
effect of surface contamination (oxide films, decomposed oils) 
on the vaporization of mercury. Some of these data are re- 
produced in Figures 4-38 and 4-39. 

It may be seen in Figure 4-38 that as the oxide layer on the 
surface is increased, the maximum temperature difference for 
the surface vaporization regime decreases . 

Many wall materials form microscopic crystalline structures 
in their surface oxides, which act as roughness elements on 
the surface. These oxide film structures are not stable on 
some refractory metals, however. If the metal bulk is not 
oxygen saturated, oxygen from surface oxides formed over a 
short period of time tends to diffuse into the metal wall and 
disappear if the temperature is high and the source of oxygen 
is removed. 
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Figure 4-35, The lifetimes for water droplets vaporizing on polished and sand- 
blasted nickel surfaces. 
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SECONDS 





Figure 4-37. The lifetimes for Freon 11 droplets vaporizing on polished and oxidized surfaces. 



LIFETIMES, 9 . SECONDS 


EX P. NO 
• 16.17 


SURFACE: TANTALUM 
LIQUID:MERCURY 
ATMOSPHERE: ARGON 
ALL OXIDATION AT 670°C 

SURFACE TREATMENT 


16.17 MECHANICALLY POLISHED. THEN CHEMICALLY 
POLISHED 

18 OXIDIZED AT 670 MM Hg IN 02 FOR I MINUTE 

19 OXIDIZED AT I00MM Hg IN OjFOR 4.5 MINUTES 

20 OXIDIZED AT I00MM Hg IN O2 FOR 5 MINUTES 

21 OXIDIZED AT 100mm Ha IN 02 FOR 6 MINUTES 

* NOT PURPOSELY OXIDIZED 

NORMAL OXIDE FILM FORMATION AT ROOM 
TEMPERATURE 


1 50 200 

T W~ T SAT'° C 


300 350 


Figure - 1 - 38 . The lifetimes of mercury droplets vaporizing on a tantalum 
surface with various amounts of surface oxidation. 
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LIFETIME, 0, SECONDS 



t w -T sat » ° c 


Figure 4-39. T..e lifetimes o' me ••.m y droplets vaporizing on a tantalum 
surface wi‘ * ■’ 'ious contaminants. 
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Table 4-1 gives an approximate range of temperature differences 

which bound the three boiling regimes discussed in the text. 

Some of the major effects of the chemical and physical features 

of the boiler tube surface are summarized below: 

1. In the surface evaporation region, wetting droplets 
evaporate as thin films and nonwetting droplets as hemi- 
spheres . 

2 . The value of the wall-fluid temperature difference at which 
the minimum in the vaporization lifetime occurs decreases 
with increasing surface roughness. 

3. The minimum droplet vaporization lifetime increases with 
surface roughness. 

4. Oscillation of droplets during film boiling increases the 
vaporization lifetime as much as 50 percent depending 
upon the type of metal surface and surface temperature. 

The height of the oscillation increases from approximately 
five to ten diameters at the beginning to as high as 500 
diameters at the end. 

5. Roughening a heat transfer surface reduces the heat 
transfer in the lower portion of the transition region 
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TABLE 4-1. 

APPROXIMATE CRITERIA FOR BOILING REGIME TRANSITION 

Smooth, clean surface Hough, clean surface 

Fluid Maximum Temperature Minimum Temperature Maximum Temperature Minimum Temperature 

difference for surface difference for film difference for surface difference for film 



vaporization 

boiling 

vaporization 

boiling 



'C 

°C 

°C 

Water 

70 

120 

55 

185 

Mercury 

70 

195 

35 

195 

Potassium 

uo 

185 

— 

— 

Freon 11 

60 

120 

35 

110 

Freon 114 

35 

55 

55 

80 


35 yo 


Benzine 



and increases the heat transfer in the upper portion of 
the transition region. 

6 . The value of the wall-fluid temperature difference at the 
maximum vaporization lifetime decreases with increasing 
surface roughness. 

7. The value of the maximum lifetime for vaporizing droplets 
tends to decrease with surface roughness. 
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DISCUSSION 


rhis chapter has presented information on transition criteria for the 
boiler design functions given in Chapter 3. This transition information 
can be used as a guide for the designer as he makes decisions on which 
phase distributions exist in particular locations proceeding throughout 
the boiler tube. Clearly, in many instances, he must use the function 
which covers the "worst case” because exact criteria on transition 
are not available 
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NOMENCLATURE FOR CHAPTER 4 



a phase distribution parameter 
tube inside diameter 
equivalent diameter of passage 
entrained droplet diameter 
droplet diameter 
helical wire insert diameter 
vapor drag force 
gravitational force 
surface tension force 


mass flow rate per unit of passage cross section 
acceleration of gravity 

proportionality constant in Newton's second law 


acceleration of gravity for earth 


normalized conductance defined by H 
heat transfer conductance 



heat transfer conductance between wall and liquid layer 
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heat transfer conductance between wall and apor 
a dimensionless parameter defined in equation set 4-23 
thermal conductivity of wall material 
axial tube length 

location in tube where "burnout? occurs 
helical boiler tube insert pitch 
heat flux 

"burnout 11 heat flux 

tube radius, helical insert centerbody radius 

liquid bulk temperature at point of first void formation 

mixed mean (bulk) temperature 

equilibrium saturation temperature 

temperature of the heat transfer surface 

liquid velocity 

superficial liquid velocity (mean velocity if liquid only flowed 
in the tube) 

superficial vapor velocity (mean velocity if vapor only flowed 
in the tube) 

vapor velocity in core of annular flow 
normalized velocity component normal to tube axis 
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x * vapor quality 

x * quality at location in tube where "burnout” occurs 

B 

z a axial distance 

* 

z = distance from initiation of heated length divided by total 

heated length 

z a axial location of first void fcrmation 

o 

GREEK SYMBOLS 


<<*> = average void fraction at an axial station 

= critical contact angle at film rupture 
6 = annular liquid film thickness 

i = tube wall thickness 


Al = a length scale defined in equation set 4-5 

AT = heat transfer temperature difference 

= defined in equation set 4-28 

4T = T - T 

s sat o 

AT , = T - T 

superheat liq ud sat 

6^ = lifetime of liquid droplet on a hot surface 
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liquid density 



P « vapor density 
v 


DIFFERENTIALS 


AM 

az 


two phase pressure gradient 
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CHAPTER 5: BOILER DESIGN EXAMPLES 


Quantitative data on forced flow boiling and pressure loss were presented 
in Chapter 3 and the corresponding flow transition criteria for the im- 
portant phase distributions were given in Chapter 4 . These heat and 
momentum transfer relations can be used by the engineer to develop 
boiler designs. 

In the following paragraphs, there are found a number of example boiler 
designs which are based on or can be analyzed with the design functions 
described in Chapters 3 and 4. The variety of applications and working 
fluids of the example boilers is one measure of the general applicability 
of the information presented in Chapters 3 and 4. A brief background for 
each boiler design is given in addition to the results obtained, including 
comparisons with actual operating performance. Detailed information 
about these heat exchangers can be found in the references that are cited. 

DESIGN EXAMPLES 

A. Design Example No. 1: Double Containment Boiler for SNAP-8 
Mercury Rankine Space Power System 

An important government program was the SNAP-8 Mercury Rankine 
cycle nuclear-electric power conversion system for space missions. 
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The power conversion system was developed by Aerojet Nuclear 
Systems Company for NASA and the nuclear reactor was developed 
for the Atomic Energy Commission by Atomics International. A 
schematic drawing of this 35 kWe system is shown in Figure 5-1. 
Thermal energy is transferred from the reactor to the power 
conversion system by a primary eu.ectic sodium-potassium alloy 
(NaK) circuit. Mercury is the working fluid 'in the Rankine cycle 
circuit. The thermal energy exchange takes place in a counterflow 
NaK-to-Hg heat exchanger where the mercury is preheated, vaporized, 
and superheated in a "once- thro ugh" boiler configuration. The 
superheated mercury vapor drives the turbine-alternator assembly 
to produce electrical power, and then is condensed and pumped back 
to the boiler. Heat given up by the mercury as it condenses is carried 
by a secondary NaK heat transfer circuit to a space radiator. Organic 
fluid in a fourth circuit is pumped to the turbine-alternator assembly 
for lubrication and cooling, and to the electrical components for cooling. 

The once-through boiler concept was selected because of its suita- 
bility for operation in zero gravity environment, and mercury was 
selected because of its modest vapor pressure at the elevated tempera- 
ture level required in order to radiate waste energy to space. Several 
versions of this boiler were built and others proposed. The design 
on which the most operating hours were obtained is described here. 

Mercury and NaK are solvents for most ordinary metals. Tantalum 
is resistant to mercury and to high purity NaK, but it was found 
impossible to maintain adequate purity level in the NaK circuit, and 
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Figure 5-1. SNAP-8 system schematic diagram. 




a bimetallic system, using stainless steel in contact with the flowing 
NaK, was required. The boiler had to be a compact, high flux design, 
and the uncertain contact resistance between layers of a bimetallic 
tube lead to the concept of two tubes, spaced apart, with the inter- 
vening space filled by stagnant purified NaK, which does not seriously 
attack tantalum. The outer tube was chosen to be oval, a shape which 
leads to the lowest average thermal resistance in the NaK layers and 
still allows for radial tube movement caused by differential thermal 
expansion when in a coiled configuration. 

in order to obtain adequate heat transfer surface and flow cross 
section, seven tubes in parallel were used and arranged in a bundle 
within a NaK containment shell. The entire boiler, which was over 
20 feet long, had to fit within a 60-inch diameter cylinder, so that 
this entire assembly was coiled. In order to make the individual 
tubes in the seven tube bundle of equal length, the tube bundle was 
twisted one full turn within the shell. The individual boiler tubes 
contained helically grooved plug inserts, followed by helical wire 
inserts. Draw ings of a few of the details of this assembly are shown 
in Figure 5-2. 

The basic method of analysis* was to subdivide the boiler length 
into a series of small mnrements and to write the coupling heat 
transfer and pressure drop relationships for each increment, 
choosing the relationships on the basis of an assumed phase distri- 
bution. The boiler preheat and superheat sections were subdivided 
by taking equal increments of vapor quality. Solution of the coupling 
equations at each node provided the local thermal and thermodynamic 
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figure 6-2. SNAP-8 bolter detail** 


conditions, and increment length. Sequential solution of the nodal 
equations resulted in an overall boiler performance prediction. 

The vapor quality-phase distribution profile, and the location of 
the inserts, is shown in Figure 5-3. Each boiler tube contained a 
36- inch long multigrooved insert plug, which extended from the 
subcooled liquid region to about 15 percent quality. This section 
was followed by a helical wire insert of p/D = 3.1, which » vt *» F i ded 
to the end of the boiler tube. Appropriate phase distributions are 
referred to in Figure 5-3, together with references to the heat 
transfer and pressure drop relationships used in these regions . 

The validity of these boiler designs was successfully supported by 

the test results of experimental single tube and full-scale tantalum- 

stainless steel boilers. Comparison of test data with design 

predictions is depicted in Figure 5-4. The analytical mercury 

temperature prediction (T^^) shows the sensible heat addition up 

to the preheat length termination point (AT ) ■ hi the two phase 

PP 

region, T^^) I s ? the vapor saturation temperature, a function of the 
local pressure, and thereafter it represents the vapor superheat 
temperature. In the boiler tests, mercury temperatures were 
measured only at the ends, so no experimental points are found 
along the curve of mercury temperature. In counterflow heat 
exchanger-type boilers the minimum, or "pinch -point, " tempera- 
ture difference lies at the point of boiling inception. 
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The first series of tantalum-stainless steel boilers (denoted by 
BRDC Nos. 1, 2, and 3), were significantly longer than the one 
described in Figures 5-2, 5-3, and 5-4. These were designed with 
a four foot long multipassage plug insert and had an overall length 
of 37 feet. BRDC Boiler No. 1 was employed in system and en- 
durance tests at a NASA facility and at a General Electric, Evendale, 
Ohio, test facility. A total test time of 15, 125 hours and six startups 
were accumulated on this boiler. BRDC Boiler No. 2 was operated 
with a breadboard power conversion system for over 8700 hours and 
accumulated 27 startups. BRDC Boiler No. 3 operation in the NASA 
test facility resulted in 135 startups and a total test time of 157 
hours. The performance characteristics of these boilers ranged 
from poor to excellent and were dependent on the boiler contamination 
state. The effects of boiler contamination ( " deconditioning" ) were 
experienced primarily in combined systems testing, and were 
attributed to oil influx into the mercury loop and decomposition in 
the boiler. This was the most serious problem encountered in the 
boiler during combined systems operation. At the operating tempera- 
ture of the SNAP-8 system, mercury wets clean tantalum and the 
decomposed organic coolant apparently caused the wall to dewet. 

A more complete discussion of the mechanisms of boiler performance 
degradation and their analytical interpretations are provided in 
Reference 2. 
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B. Design Example No. 2: Boiler for Potassium Rankine Space Power System 

NASA Lewis Laboratory has had significant interest in space Rankine 
cycle power systems where the working fluid is potassium. As a 
result of this interest, a test facility for potassium Rankine cycle 
power plant development was constructed and operated by the Space 
Power Division of the General Electric Company.** The system 
consisted of three flow circuits interconnected by heat exchangers. 

These were lithium-heating, potassium-boiling, and NaK heat- 
rejection circuits. The lithium flowed from an electromagnetic 
pump through a flowmeter, a 400-kW capacity AC electrical heater 
which simulated a nuclear reactor, and into the shell-side of the 
single-tube boiler; from the exit of the boiler, the lithium was returned 
to the suction side of the pump for recirculation. Liquid potassium 
flowed from a second electromagnetic pump through a flowmeter, a 
liquid throttle valve, a 30-kW capacity AC electrical preheater, and 
into the tube side cf the boiler. From the boiler, the potassium 
vapor entered a radiant desuperheater and vapor throttle valve which, 
in combination, simulated the turbine of the power plant. Expanding 
through the throttle valve, the vapor then flowed into a three-tube 
condenser, where heat was removed from the vapor and transferred 
to the condenser coolant, NaK. The potassium in the liquid state then 
returned to the potassium pump for recirculation to the boiler. The 
NaK coolant flowed from a third electromagnetic pump through a flow- 
meter and into the shell-side of the condenser. From the condenser, 
the NaK entered an air-cooled heat exchanger and then returned to 
the pump. 
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Both the lithium and potassium circuits were fabricated of the 
tantalum-based refractory alloy, T-lll, and were contained in a 
vacuum chamber about 5.5 meters tall and 1.8 meters in diameter. 
This chamber was capable of extremely low pressures and was 
equipped with bake-out heaters and cooling coils. The NaK loop 
was fabricated of Type 321 stainless steel. Most of this circuit 
was outside of the vacuum containment. 

A schematic diagram of the boiler is shown in Figure 5-5. The 
curved test boiler consisted of an outer and inner tube, both ap- 
proximately 2.29 meters long, to which were attached appropriate 
fittings for welding into the liquid metal circuits. The inner tube, 
in which the potassium vaporized, had an inside diameter of 1.7 
cm. The lithium flowed countercurrent to the potassium in the 
0.45 cm annulus formed by the two tubes. Pressure gauges 
and thermocouple immersion wells were located near the entrances 
and exits of the lithium and potassium passages. The boiler tube 
contained a composite swirl-generating insert for enhancing the 
heat transfer to the potassium. The initial 1.4 meters of the insert 
consisted of a single helical vane wrapped about and spotwelded to 
a 0-635-cm-diameter hollow centerbody. The remaining 0.89- 
meter length of the insert consisted of a helical wire coil. The 
ratio of pitch to the tube-inner-diameter for the composite insert 
was about 3.0. The centerbody of the helical-vane portion of the 
insert served as an immersion well into which seven thermocouples 
were placed for measuring potassium fluid temperatures. Seven 
additional thermocouples were placed in a similar tube in the region 
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Figure 5-5. Schematic drawing of single-tube boiler. 



of the boiler which contained the wire coil. The lithium tempera- 
ture distribution was measured by four thermocouples distributed 
uniformly around the shell periphery at each of 15 axial stations . 

Figure 5-6 presents experimental temperature data for this boiler, 
for the inlet and outlet conditions given in Table 5-1. Table 5-1 also 
lists the inlet and outlet conditions for the boiler at the nominal 
design point of the advanced Rankine system. From this table, it 
is clear that the run shown in Figure 5-6 closely matches the nominal 
operating condition required for the system. From the temperature 
profile shown in Figure 5-6, the liquid, two-phase, and superheat 
regions of the potassium can readily be identified. The liquid region 
extends from the beginning of the heated zone to a point about 0.15- 
meter downstream of the inlet. Within this very short length, the 
potassium temperature rises from about 1040 °K to the local satu- 
ration temperature of 1387 °K. The two-phase region starts at the 

0. 15-meter point within the heated zone and extends to the 1.5-meter 
point. From the end of this region to the boiler exit is the super- 
heated-vapor region. 

Several phase distribution regimes would be expected in this boiler. 
These are: single phase liquid superheating, wetted wall annular 
flow, and entrained flow with a dry wall. The latter phase distri- 
bution probably extends into the net vapor superheat region where 
superheated vapor carrying droplets of saturated liquid exist in such 
proportion that the mixed mean temperature is superheated. The 
design relationships for flow through a helical vane and helical wire 
would be appropriate. In the high quality and superheat regions. 
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Figure 5-6. Boiler temperature distribution. 
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TABLE 5-1. 


Comparison of Boiler Inlet and Outlet Conditions 
At the Nominal Operating Point of the Potassium 

3 

Rank me System With Those for Run Shown in Figure 5-6. 


Parameter 

Nominal Operating 
Point of 

Potassium System 

Test Conditions 
for Run Shown 
in Figure 5-6. 

Lithium inlet temperature, X 

1477 

1456 

Potassium inlet temperature, K 

928 

1039 

Potassium exit temperature, K 

1422 

1*23 

2 

Potassium exit pressure, N/cm 

112 

109 

Potassium exit superheat, K 

40 

47.8 
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the functions for helical nonwetting droplet flow would be 
applicable. 

A plot of the pressure losses across the boiler as a function of 
potassium flow rate is given in Figure 5-7 . These data were 
obtained during runs in which the potassium-exit pressure (exit 
saturation temperature) and the potassium-inlet temperature were 
held constant. The flow rate and the heat input were varied to 
achieve an exit quality of 1.0 and two levels of exit superheat. 

The pressure losses, as shown in this figure, are substantially 
larger for the boiler operating with a modest exit superheat. This 
increased pressure loss is due to the longer vapor path lengths 
present within the boiler tube and to large momentum pressure 
losses. 

C. Design Example No. 3: Rankine Cycle Boiler for an Organic 
Working Fluid 

The Army has undertaken the development of several Rankine 

cycle engines utilizing organic working fluids. These engines are 

intended as silent prime movers for electrical power generators. 

4 

The boiler for one of these engines is described in the following 
paragraphs. 

Because of the low noise requirement, the combustion heat source 
for this boiler was of unusual design and recycled a portion of the 
exhaust gases into the intake. The burner was, therefore, integrated 
with the boiler. A typical burner/boiler module is shown in Figure 5-8. 
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Figure 5-7. Effect of exit quality and vapor superheat on 
boiler overall pressure loss. 
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In addition to the burner and primary heat exchanger, the unit also 
contained control equipment and instrumentation appropriate to this 
application. Figure 5-9 is a schematic diagram of the module. 

The organic vapor temperature was monitored by a thermocouple 
p».#be. which operated a temperature controller. The temperature 
controller regulated the combustion rate (by means of the fuel pump 
motor speed) to maintain a constant organic vapor outlet temperature 
of the boiler, despite variations in organic fluid flow rate. 

The heat exchanger consisted of two staggered tube arrays of 
stainless -steel finned stainless steel tubing. The lower of the 
two banks, next to the combustor, was the boiler, and the preheater 
was directly above it. The tube size was somewhat larger in the 
boiler unit to accommodate vapor flow without excessive pressure 
drop, and a single-tube design was used. The preheater section 
was a counterflow unit and the boiler was a parallel flow unit. 

Parallel flow was found to be desirable in limiting tube wall tempera - 
atures at the boiler exit. A drawing of this array is shown in Figure 
5-10, and a photograph in Figure 5-11. The important dimensions, 
calculated heat transfer rates, conductances, and pressure drops 
are given in Table 5-II. 

The organic flu . used in this unit, mono-isopropyl-biphenyl, is 
temperature sensitive, and it was determined that if the tube wall 
temperature was limited to 850° F, adequate system lifetime was 
attained. It may be seen from information in the table that the 
calculated highest temperature was on the front side of the final 
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Figure 5-9. Burner/boller functional schematic. 





COMBUSTION GASES 


Figure 5-10. Section through organic fluid preheater 
vaporizer assembly. 
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TABLE 5 -II. 


CHARACTERISTICS OF ROCKETDYNE ORGANIC FLl ID 
PREHE ATER- BOILER- FINA L CON FIGURATION 


Characteristics based on gas flow rate 13G lb/ hr 

Gas temperature 2000" F 

Organic flow rate 290 lb hr 

Organic saturation temperature TOO^F 
Organic outlet condition 10 F superheat 


Thermal Characteristics 

Boiler 

Preheater 

I. Temperatures >°F) 



Organic inlet 

700 

568 

Organic outlet 

Maximum wall temperature, organic 

710 

700 

side 

7 s7 mean outlet 
831 outlet front side 

— 

Outlet gas 

1245 

700 

2. Temperature Different s (’F) 



Gas to organic, mean 

920 

340 

Gas to organic, maximum 

1300 

545 

9 

3. Heat Fluxes (Btu hr ft”) 



Mean 

33,300 

14,700 

Maximum 

45,100 

— 

Minimum 

23,200 

— 

2 

-1 Conductances (Btu nr ft” F) 



(Eased on inside area except as noted) 
Inside tube wall t> organic, mean 

1200 

313 

Inside tube wall to organic, minimum 

266 outlet 

286 

Gas to organic, mean 

37 . 9 

56.7 

Gas t< organic minimum 

34.1 

55. 8 

Gas side to fins, mean (based on gas 



side area) 

>.s3 

15.68 
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TABLE 5-n (continued) 



Boiler 

Preheater 

B. Flow Characteristics 
1 . Pressure Drops 

Gas side 

0.025 in/H 2 0 

0.112 in/HgO 

Organic side 

6.5 lb/ in 2 

0.84 lb/in 2 

Mini mu in gas side flow area 

0.128 ft 2 

0.087 ft 2 

C- Dimensional Characteristics (inches except 
where noted) 

Tubing O. D. 

1/2 

3/8 

Tubing wall thickness 

0.020 

0.016 

Fin diameter 

0.900 

0.650 

Fin thickness 

0.018 

0.018 

Fin pitch (fins per inch) 

6.5 

10.5 

o 

Heated area (gas side) 

4.46 ft 

7.20 ft 

Heated area (organic side) 

0.907 ft 2 

1.38 ft 2 

Heated length 

7.47 ft 

15.4 ft 

Total length 

11.2 ft 

23.3 ft 

Number of tubes (total) 

16 

33 

Number of rows 

4 

5 

Arrangement 

Staggered 
4-4 -4 -4 

Staggered 
7 -6-7 -6-7 


Weight of tubing and fins 10.46 lbs 

Weight of tube sheets and side 

panels 3 . 60 lbs 

Weight of housing and flanges 4.024 lbs 

Total weight 18 -2 lbs 

Material: Stainless steel. All flow passages type 347 with nicrobrazed joints. 
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row of tubes in the boiler. The thin walled stainless steel tubing 
did not conduct a significant amount of heat around the periphery 
of the tube, so that nonuniform heat transfer from the gas to the 
fin was not equalized. This local hot zone was determined to be 
acceptable for the application, which required only a relatively 
short service life. Conductances were computed for this preheater- 
boiler system on the basis of turbulent single phase flow for the 
preheater, wetting annular flow with a turbulent film in the two 
phase region from boiling inception to 85 percent quality, single 
phase flow vapor heating with a variable flow rate to 100 percent 
quality, and single phase flow of the vapor in the very short super- 
heat region. 

The computed local two phase pressure gradient is shown in 
Figure 5-12, together with the gradient as predicted by the 
Martinelli correlation. 

Operating experience with this unit indicated that its performance 

4 

was quite close to the design calculations. 
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Figure 5-12. Two-phase flow pressure gradient in organic boiler. 
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D. Design Example No. 4: A Liquid Hydrogen Boiler 

A small liquid hydrogen vaporizer using a flow of liquid hydrogen 
as the heat source was designed and constructed for a special 
application. The arrangement employed was a counterflow heat 
exchanger, with vaporization taking place in a central tube sur- 
rounded by an annular passage containing the hot side fluid. The 
general requirements for this boiler are given in Table 5-m. 


TABLE 5- ID. 

General Characteristics of Hydrogen Boiler 



Cold Side 

Hot Side 

Flow rate 

0.44 gm/sec 

150 gm/ sec 

Inlet temperature 

17.2 °K 

20.8 

Outlet temperature 

20°K 

20.7 

Inlet pressure 

0.34 atm 

1 . 16 atm 

Pressure drop 

0.027 

0.0027 atm 


239 



The heat flows required to accomplish this performance are 
200 watts in the boiling section, and 13 watts in the superheater. 

It may be shown by means of standard heat exchanger calculations,* 
that the specified conditions require an overall conductance-heat 
transfer area product for the vaporization section of UA o 47.1 
watts/ °C, and for the superheat section of UA =6.0 watts/°C. 

Liquid hydrogen is a low surface tension fluid which wets metal 
surfaces readily. For the small overall temperature differences 
in this application, no possibility of film boiling existed. The 
phase distributions expected were turbulent liquid preheating, 
annular flow throughout most of the vaporization region, and dry 
vapor in the superheater. Because of the very small vapor ve- 
locities and heat flows, no entrainment was expected, so that the 
wetted wall condition would extend all the way to 100 percent 
quality. Round 0.81 cm inside diameter tubing was chosen for the 
boiler tube. 

The hot side duct was an annular passage with an equivalent diameter 
of 2.5 cm. In Table 5- IV are presented the principal calculated 
performance data for this unit. 


* 

See, for example. Reference 5. 
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TABLE 5- IV. 


Computed Performance of Hydrogen Vaporizer 


Vaporization in annular flow at 70 percent quality 
Void fraction 0.87 

Liquid film thickness 0.25 min 

Heat transfer conductance: 

Based on turbulent film, no nucleation 
Based on laminar film, no nucleation 


Liquid preheating: 

Heat transfer conductance 
(turbulent flow. Re = 4400) 

Superheating of vapor: 

Outlet vapor conductance 
(turbulent flow. Re = 54,000) 

Outlet vapor velocity 

Integrated two phase pressure drop, 
for annular flow 

Hot side conductance (turbulent 
flow. Re = 455,000) 


23.9 watts/cm 
2 

3.7 watts/cm 

2 

2 . 6 watts/cm 

2 

2 . 0 watts/ cm 
2380 cm/sec 

0.012 atm 

2 

31.1watts/cm °C 



It is rather difficult to test a heat exchanger under (he conditions 
specified for this design. The temperature differences and flows 
are very small, and accurate measurements of liquid hydrogen flows 
are not simply made. Therefore, a proof test of this unit was made 
using Freon 11 (another wetting fluid) as the vaporizing flow, and 
warm water as the hot side fluid. It was possible to perform these 
heat exchanger tests with very large water flows, so that the mean 
vaporization side conductances could be extracted from the data. 
Results of these tests with a Freon 11 flow and temperature differences 

chosen to match inlet Reynolds number and heat flow, yielded a mean 

/ 2 

vaporizing conductance of 29.6 watts/cm °C. The predicted mean 

vaporizing conductance for Freon, utilizing the model of a turbulent 

2 

annular film without nucleation, w *s 26.1 watts/cm °C. 
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E. Design Example No. 5: Rankine Cycle Automobile Engine Vapor 
Generator 

The design techniques discussed in Chapters 3 and 4 are also being 
used to design vapor generators for "steam" automobiles now under 
development. The working fluids for these engines include several 
very specialized organic fluids and organic fluid mixtures, and 
water. 

This section is divided into two parts, which are; 

1. A general discussion of automobile vapor generator design 
considerations, and 

2. A specific example. 

1. Rankine Cycle Automobile Vapor Generator Considerations 
Some of the considerations which the vapor generator designer 
working on this application must take into account are itemized 
below: 

(a) The timelag from an increase in water flow to an increase in 
steam flow must be short, in order that the vehicle can respond 
rapidly to changes in power demand. There must be no signifi- 
cant energy storage in a boiler for this application so that 
water flow, steam flow, and heat input (combustor output 
level) must match closely, even during rapid transients. 

Some modern designs use several parallel paths in the two 
phase and superheater portions of the vapor generator to 
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minimize transit time in this region. The multitube design 
shortens response time but can introduce .'low instabilities. 
The Doble steam automobiles, which were manufactured 
commercially for many years, used a monoti.be boiler 
with a great length of tubing and a long flow transit time- 
These boilers employed a "normalizer" system to more 
closely match steam conditions with steam demand than 
would have been p issible with water input to the preheater 
only. The normalizer was a second water input near the 
inlet to the superheater which introduced water to trim 
superheater output conditions. This water flow, being 
much closer in both distance and time to the expander 
inlet, provided the necessary fast response. This system 
is also used in the steam powered bus constructed by 
Brobeck Associates for the State of California. 

(b) The heat fluxes in the vaporizing and superheating sections 
must be chosen so that tube wall temperatures are well 
controlled. For example, it has been shown in Chapter 4 
that a large change in boiling side heat transfer conductance 
occurs when the wall is no longer wetted by the working 
fluid, and the vapor quality at which this occurs depends 
upon a number of complex interrelationships ir. geometry, 
flow rate, and fluid properties. Although an automotive 
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steam generator must be a high average heat flux 
design for compactness, the transition cited above and all 
the regions downstream must be in a low flux region to 
prevent tube failure. 

(c) Stored energy and liquid holdup in the boiler must be small 
to prevent serious accidents in the case of tube rupture or 
vehicle co'lision, and also to minimize start up time. Fuel 
economy on short runs can be seriously affected by energy 
required to bring the boiler to operating temperature. 

(d) Lubricants finding their way into the boiler must not be 
decomposed on the heat transfer walls to form thermally 
insulating or nonweuing layers. A typical historical 
solution to this problem has been to limit wall temperatures 
to levels which the lubricants can tolerate. Other solutions 
proposed to keep oils from the boiler include oil separators 
at the boiler inlet, and dry lubricants, or self-lubricating 
surfaces in the expander and pump. 

2. A Specific Automotive Vapor Generator Example 

A Rankine cycle automotive engine vapor generator for an 
engine of approximately 100 bhp is described in the follow- 
ing paragraphs. 
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As was commented on in item (b), heat fluxes must be con- 
trolled to prevent tube wall overheating. In many modern 
designs, the technique for accomplishing this is to place the 
vaporizer closest to the combustor, with the low quality 
vaporization region in cross -parallel flow, then introduce the 
suoerheater into the hot gas stream, operated cross counter- 
flow, follow with a few rows of high quality vaporizer section 
(sometimes called the "dryout region"), and place the liquid 
preheater farthest from the combustor. Such a flow path is 
shown in Figure 5-13. The example design presented here 
approaches the problem of heat flux control in a different 
manner, specifically, by employing tubing coating where 
necessary with a thermal insulator. Thermally resistive 
coatings have been utilized in portions of organic boilers 
presently in service to protect the working fluid, so the general 
technique has been tested successfully. The adoption of such 
a system for the example water vaporizer allows the use of a 
simple counterflow design. The flow path and typical conditions 
are shown in Figure 5-14. Table 5-V presents some salient 
characteristics of this design. 

It may be seen that with a superheater of two channels, that the 
maximum wall temperature is in the order of 632 °C. If, for 
strength and corrosion resistance, the outlet wall temperature 
were to be limited to 593 , 'C (for example), some thermal insulation 
would be required at the superheater outlet. This could take the 
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Figure 5-13. Flow path in typical modern automotive forced convection 
vapor generator. 
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Figure 5-14. Flow pHth and stream conditions in example boiler. 



TABLE 5-V 


Combustor heat rate 

0. 525 x 10 6 

watts 


Steam heat rate 

6 

0.42 x 10 watts 


Thermal efficiency 

80 percent 



Characteristic 

Superheater 

Vaporizer 

Preheater 

Heat transfer effectiveness 

0 208 

0.581 

0.518 

Arrangement 

multipass 

cross 

counterflow 

multipass 

cross 

counterflow 

multipass 

cross 

counterflow 

Stream capacity rate ratio 
C 

min 




C 

max 

0.75 

0 

0.62 

Conductance-heat transfer area 
product UA watts/ °C 

51.6 

280 

274 

Gas side conductance 
watts /cm 2 *C 

0 . 0169 

0.0159 

0.0159 

Water side conductance 
watts/cm 2 °C 

0.133 

variable 

with 

location 

1.42 

Overall conductance 
(without insulation) 
watts/cm 2 °C 

0.0153 

0.0142 

0.0156 

Maximum temperature 
difference °C 

816 

843 

343 

2 

Heat flux watts/cm 

16.4 

15.6 

7.2 

Wall temperature C 

632 

393 

289 

Tube diameter (bare) mm 

19.1 

19. 1 

12.7 

Tube spacing transverse mm 

23.9 

23-9 

16.0 

Tube spacing longitudinal mm 

23.9 

23.9 

16.0 

Number of parallel channels 

2 

1 

1 
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form of a ceramic coating, or a metal foil wrap spaced from 
the tube by dimpling, or some other simple device. Inspection 
of the wall temperatures row by row would be necessary for 
this arrangement, particularly under off design or transient 
conditions. 

This design, utilizing a single vaporization path and two 
parallel superheater paths, is not particularly compact. A 
significant saving in tube length can be realized in the cooler 
portions of the assembly by use of finned tubing, rather than 
the close-spaced bare tubing of the example. 

This application is interesting in that the gas side conductance 
controls the heat flux and the boiling conductance variations do 
not affect the design. Only the superheated vapor conductance 
needs to be computed accurately. This is not the case for the 
geometry in which the vaporization takes place in the tubes 
next to the combustor. For that case, boiling conductance 
heat fluxes and wall temperatures must be evaluated carefully 
to prevent tube overheating. 
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CONCLUDING REMARKS 


It has been shown that the heat transfer and pressure drop functions 
presented in Chapters 3 and 4 have been used successfully to design 
boilers whose measured performance characteristics are in satisfactory 
agreement with predicted values. 

Results of the NASA SNAP-8 mercury boiler program have shown that if 
the mercury is kept free of foreign materials such as oxygen and organics 
and if optimized designs (based on helical flow) are used, mercury boilers 
with good heat transfer and fluid flow performance characteristics can be 
built and successfully used in industry. It has also been demonstrated in 
NASA's potassium Rankine cycle program that efficient potassium boilers 
can be designed, fabricated and operated. In addition, working fluid purity 
levels required have been established. These high temperature boilers can 
be used to extend the thermodynamic efficiency beyond that attainable with 
mercury. 

The major limitation in organic boilers is the degradation of these fluids 
at elevated temperatures. When an organic working Quid is exposed to 
high temperature levels, the chemical composition of the Quid changes and 
the degree depends upon container metal, time, as well as temperature 
levels. The boiler designer must consider these processes and the system 
thermal efficiency to determine how they inOuence the overall economics 
involved . 
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The major reason for the development of a Rankine cycle automobile engine 
is that the external combustion process can be optimized to reduce generation 
of the pollutants formed in an internal combustion engine. The Environmental 
Protection Agency, the State of California, and several industrial organizations 
are constructing Rankine cycle automotive engines using a variety of working 
fluids. 
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CHAPTER 6: BOILER APPLICATIONS 


Five example boiler designs were presented in Chapter 5 to illustrate 
how the design functions and transition criteria presented in Chapters 
3 and 4 can be utilized. These examples also illustrate some appli- 
cations of forced flow and once-through boilers to technologies other 
than space technology. It is possible to envision additional boiler 
requirements and applications, some of which are summarized in 
the following paragraphs. 

A . Mobile Propulsion Systems 

An example of a forced flow boiler application to future mobile 
propulsion systems (namely, the automotive steam car) was noted 
in the previous chapter. Interest has been developed in this system 
because -he external combustion process associated with it appears 
to have excellent potential for reducing combustion gas emissions. 

It is believed that mobile power plants which utilize forced flow 
boilers could also be used in trains, trucks, and buses. In such 
systems, high specific power densities and good operating per- 
formance characteristics can be achieved. 

B. Topping Cycles for Power Plants 

The concept of using a topping cycle in a power plant to increase 
the overall thermodynamic efficiency is an old one. However, this 
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approach is now more appealing because the technology of high 
temperature, low pressure working fluids (such as tue liquid 
metals) has now been developed in addition to the forced flow and 
once through boiler concept. Potassium i3 of particular interest 
because *t has high temperature capabilities, has good thermal 
properties, and its cost is not excessive. 

C. Bottoming Cycles for Power Plants 

It is thought that there may be some advantages in including a 
low temperature Rankine cycle system in a power plant to utilize 
some of the associated waste heat. There are a number of organic 
working fluids whose thermal properties are sufficiently good and 
whose decomposition rates are negligible at low temperatures so 
that such bottoming cycles may modestly increase overall power 
plant efficiencies. 

D. Auxiliary Power Plants 

As new land areas throughout the world are developed (particularly 
in remote regions where it is difficult and expensive to route electric 
power lines), it will probably be important to construct auxiliary 
power plants whose heat source may be either nuclear or fossil 
fuel types. In order to keep the capital costs for such systems 
to a minimum, the designers will clearly desire to have compact 
systems; optimized boilers would be of particular importance. 
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E. Advanced Refrigeration System 

In order to r- duce the fuel and energy requirements for large 
refrigeration systems, more advanced refrigeration systems will 
be required. Most modern refrigeration system evaporators 
operate on a semi-forced flow basis in that the boiling process 
is more accurately defined by a natural convection type re- 
quiring large volumes. It is believed that more efficient and 
compact evaporators, based on forced flow boiling, may help in 
obtaining increased system efficiency and reduced component 
sizes. 

F. Boiling Nuclear Reactors 

Nuclear reactor designers have considered the advantages of 
systems that contain boiling working fluids. Boiling water 
reactors have, of course, already been built and are in operation. 
It is felt that the use of forced helical flow boiling in such nuclear 
reactor cores will not only yield high power density systems but 
reduce the possibilities of slug flow thereby increasing fluid flow 
and reactor stability. 

G. Solar Evaporation Boilers 

There is a great interest in utilizing solar energy as a prime 
source because of the national objective of conserving fossil 
fuels. One type of conversion process involves the evaporation 
of a working fluid, which is then used in a Rankine cycle system. 
In such a concept, forced flow boiling again would have a signifi- 
cant advantage. 
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H. Geothermal Power Systems 

The government has identified our geothermal heat sources as having 
great energy value. Hot fluids pumped through fractured geothermal 
strata would change phase at points of pressure reduction near the 
surface of the earth where wet steam would be transferred to some 
type of expander or turbine. In order to prevent mechanical erosion 
by the wet steam, a high vapor quality boiler and superheater would 
have to be interposed between the geothermal steam source and the 
prime mover. Again, forced flow boiling in this component would 
be used. 

I. Process Heat Exchangers Involving Change of Phase 

The process industry utilizes a variety of evaporative heat exchangers 
and many different working fluids. The forced convection boiling 
technology presented here can be applied to such heat exchangers 
to provide better temperature control (by preventing hot spot regn.s) 
and to reduce capital costs (by optimizing size and pumping power 
requirements) . 
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APPENDIX A: THERMAL PROPERTIES OF SOME OF THE 


IMPORTANT WORKING FLUIDS 


Figures A-I through A -7 contain thermal conductivity, specific heat, 
viscosity, density, Praudtl number, vapor pressure and surface tension 
data as functions of temperr ture level for some of the important working 
fluids in the liquid and saturated vapor states. These data are presented 
for the convenience of the boiler designer. 
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Figure A-5. 
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VAPOR PRESSURE, KG/ CM 



Figure A-6. 
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Figure A-7 
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APPENDIX B: 

EQUATION SETS FOR CHAPTER 3 

The experimental semi-empirical and analytical information presented 
in tiie figures in Chapter 3 are defined by functions that represent the 
heat transfer and pressure drop design functions in this monograph. 

These functions are presented in this Appendix in the form of equation 
sets that are numbered in conformance with the figures in Chapter 3. 

Equation Set 3-1: 

64 

£ = — , analytical and experimental 

0.023 

£ = — , experimental 

Re 

Equation Set 3-2: 

Turbulent single phase 

0.8 0.4 

Nu = 0.023 Re Pr , experimental and semi-empirical 
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Equation Set 3-3: 


Turbulent single phase 

0.8 

Nu = 7 + 0.025 Pe , analytical 

0.4 

Nu = 0.625 Pr , experimental and semi-empirical 

Equation Set 3-4: 

Turbulent single phase (analytical) 


Nu 


F 

o 


rr at 


Vf 


c Re Pr 


8 AT 


m 


5 |y Pr + In (1 + 5yPr) + 0.5 F q In J 


A function of Y, Pr, Re, and C • 


Equation Set 3 -6: 


Nonwetting slug flow 


^ P TP 
A p (x = 1) 


exit 


(f- 


2 x 


, analytical 


exit 


26 8 



Equation Set 3-8: 


Nonwetting slug flow 
Laminar flow 


H 




, analytical 


Turbulent flow 



0.8 

+ x 


1 



V 


X 


0 


X 


, analytical 
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Equation Set 3-9: 


Wetting slug flow 


a 


u 

v 


1.2 (u. + u ) + 0.35 
I v 


VgD 


, analytical 


Equation Set 3-10: 


Upward annular flow , analytical 



r u 

8 It / o v 
t \ w 

V 


)(; 


w u 
I 1 


w tl 


V V 



5 


1 

4 



o 




o 
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Equation Set 3-14: 


Viscous-turbulent annular flow, analytical 
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Equation Set 3-15: 


Void fraction versus pressure drop 


itt = T~^\t rz^7s * analytical * Reffcrence 19+ 
v + r yF) 


<t> 

i « 



analytical, Reference 22* 


Equation Set 3-16: 


Turbulent annular film 


c 

w D 
I 

k = 0.0131 

t nr 


1 + 1 


•*e) 


, analytical 


★ 

References are listed at the end of Chapter 3. 
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Equation Set 3-17: 


Wetted rivulet flow, analytical 


h 

v 
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Equation Set 3-19: 


Annular flow with entrainment, analytical 



o 


( 4 ) w f = f* t F(y + ) 

o 
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<5) 


F(y + ) = 3y + + 2.5y + In y + - 64 for y+ > 30 

F(y + ) * 12.5 - 8.05y + + 5y + In y + for 5 < y + < 30 

F(y + ) as 0.5y +2 fory + < 5 

> 
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Equation Set 3-23: 


Single phase turbulent 




, semi-empirical 


Equation Set 3-24: 

Helical flow single phase 




— , analytical 


Equation Set 3-25: 

Helical flow single phase 


Nu 


forced 


Nu. 



, analytical 
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Nu 


free 


Nu 


Re 013 Pr 0 #7 ' p/D ' ' D ) ' 


Nu. 


Nu forced "“free 
Nu. + Nu 


analytical 


D = D = D 
I h i 


Equation Set 3-26: 

Helical wetted rivulet flow, analytical 


h 

vap 



(1 - X) 


8/15 


x 


1/5 



0.4 , 0.6 1.2 

6.6 L k w 

v 


p 02 p 021 

I V 




(g) 


0.4 
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Equation Set 3-28: analytical: 



m = 1 for twisted tape 
m = 2 for helical vane 
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Equation Set 3-32: 


Helical wires single phase, analytical. Reference 25 



Foi round wires, ~ 1.0 

w 
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Equr'ion * >t 3-33: 

Hjlical wires single phase, analytical, Reference 25 

v 
v 

without slip 

K is defined in Equation Set 3-32. 
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APPENDIX B: 


EQUATION SETS FOR CHAPTER 4 


The experimental, semi-empirical, and analytical equation sets for 
Chapter 4 are presented below. 

Equation Set 4-3: Analytical 


+ 



D. 




f P. 


2 

u D 
I d 


32 g ^a<t> 


D d( p < : »v) e 

12 <J<f> g 
c 


For y , <5 
d 


T 

sat 



q 


c p 
P f 




5 In (1 + 5 Pr) 


+ 
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For 30 < y + 
J d 


T - T 

sat B ! l 2.5 

[ 1 
+ 

y d ! 


L'i 1 

q l°P P < “«l = JI 

8 

+ 

y CL 


y CL 


0.02 <<t> < 0.03 


Equation Set 4-5: 


Semi-Empirical 



Length required to heat liquid to saturation in the absence of boiling 


At = i P. u. c (T - T ) 

f q I I p s o 
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Equation Set 4-6: 


«*>= 


P, L A 
I c 


(z - z ) - u 
o I . 


0 . c 4T 
I p s 

P L 
v 


) q /*h \ l i P t C d 4T s MK-'I i H c 4T (1 - F )i 

P~L A~ < z - 7 >- P — V 1 * C u ♦ F ''l + -* S. L 

< P I c o' P» v L )( p J \%S i F 2j{ ,+ L j 


F 


1 



(from Equation Set 4-5) 


F 

2 


1.18 


ag (P ! " P v ) 


1/2 


The initial point of vapor generation for water has been found •mpirically 
to be approximated by 


T - T 
s o 

q 



14 + O.lp 


cm 

where u is given in 

*■ 4 sec 


and p is »n atmospheres. 
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Equation Set 4-7 : 


Semi- Empirical 




* 

u 

V 


(7 + 0.06-^- ) 



P 

v 


u*<10 

1 


* * 

u = 0.9 + 0.6 u. 

v l 


Equation Set 4-11 


Analytical 


V 

v 



> t 

g e P °p + P D 

f V V D 

3 

C D 

V 


1/2 
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Equation Set 4-18: 


Analytical 


(a)*" 


Equation Set 4-20: 


Analytical 



a - - sin 2 a 
c 2 c 

2 

(sin a ) 
c 


Equation Set 4-21: Analytical 

(1) Viscous -viscous annular flow (Chapter 3) 

(2) Equation Set 4-20 
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Equation Set 4-23: 


Analytical 


v 

t 
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Equation Set 4-26: 


Analytical 


Small Slip Model: 


2 

_ -E- 

[ 2 c - 2 - 

h. 

1/2 r .2 r > 1/2 

1 -( ° 

D 

D 7T 

8 D D 

p K 

1 \ D / 


P 

£ 

1- 


The iarge slip model is not easily represented in closed form. See 
Reference 5 for computational details. 


Equation Set 4-28. 

The analysis does not yield a result in closed form. The definition of 
AT l is 


"i ■ 


( T i ) - ( T , 


i^ide wall 


at -x) 


outsiue wall 


H T i 


inside wall 


)l 
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Equation Set 4-34 


Analytical 


Surface vaporization (constant contact area). Reference 35. 


0 . 


8 < i Lr » 6 c /V \ 4 

9 k | <*« - V f c « ' r d' 


Surface temperature for minimum vaporization lifetime for surface 
vaporization 


2 o cos a 
AP V = r 
P 


Change in lifetime with spallation 


0 , 




n = 1 


( volume of divided droplet 
volume of undivided droplet 


7/15 
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Film boiling on a smooth surface 


e. = 


30f 


1/5 


4/5 _ 3/5 7/5 

'l L 'd 


1/5 1/5 „ 

7(160) p (T 
v w 


T 3/5 k 3/5 1/5 

S V 


The relationships for oscillating droplets and for film boiling are not in 
closed form. Derivations and results may be found in Reference 35. 
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NOMENCLATURE FOR APPENDIX B 


flow cross section 
a phase distribution parameter 
drag coefficient 

drag coefficient for a wire against a wail, modified for 
tube boundary layer 

liquid specific heat 

tube inside diameter 

vapor bubble diameter at departure from wall 

equivalent diameter of passage 

entrained droplet diameter 

droplet diameter 

helical wire insert diameter 

functions defined in equation set 4-6 

functions defined in equation set 4-6 

vapor drag force 

gravitational force 
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surface tension force 


friction factor 

fraction of droplet base area in direct contact with surface 
mass flow rate per unit of passage cross section 
acceleration of gravity 

proportionality constant in Newton's second law 
acceleration of gravity for earth 

h 

normalized conductance defined by H = 7 — r 

h = ioOtd 


x 


heat transfer conductance 

heat transfer conductance between wall and liquid layer 
heat transfer conductance between wall and liquid vapor 
a dimensionless parameter defined in equation set 4-23 
thermal conductivity of liquid 




thermal conductivity of vapor 
thermal conductivity of wall material 
latent heat of vaporization 
axial tube length 

normalized length defined in equation set 4-23 
axial length in tube 

location in tube where "burnout" occurs 
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> in 



Prandtl modulus of liquid 
helical boiler tube insert pitch 
heat flux 

" burnout" heat flux 

tube radius, helical insert centerbody radius 
radius flat disk with volume equal to that of round drop 
initial drop radius 

meniscus radius of liquid-vapor interface in surface cavity 
radius of round drop 

liquid bulk temperature at point of first void formation 

temperature on the heat transfer wall at a conductance 
discontinuity 

mixed mean (bulk) temperature 
equilibrium saturation temperature 
temperature of the heat transfer surface 
liquid velocity 
liquid inlet velocity 

superficial liquid velocity (mean velocity if liquid only flowed 
in the tube) 
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superficial vapor velocity (mean velocity if vapor only flowed 
in the tube) 


* 

u 

v 


V 

v 



W 

t 

x 


X 


B 



z 


vapor velocity in core of annular flow 

transverse velocity component at the wall in helical flow, 
assuming solid body rotation 

normalized velocity component normal to tube axis 

helical flow transverse velocity component for fully established 
flow 

total mass flow in tube 
vapor quality 

quality at location in tube where •burnout" occurs 

y + evaluated with the tube radius rather than the bubble 
radius 

defined in equation set 4-3 
axial distance 


GREEK SYMBOLS 
a = contact angle 

<a> = average void fraction at an axial station 
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a = critical contact angle at film rupture 

c 

6 = annular liquid film thickness 

6 = tube wall thickness 

Al = a length scale defined in equation set 4-5 

/ip = excess pressure above liquid saturation pressure at local 
temperature in cavity 

AT = heat transfer temperature difference 

A^ = defined in equation set 4-28 

AT = T - T 

s sat o 

y^fp , fj 1 ^ rp 

superheat liquid sat 
£ = friction factor for a smooth tube 

6^ = lifetime of liquid droplet on a hot surface 

Ol ^ = lifetime of one of n equal droplets resulting from spallation 

k = a parameter in equation set 4-23 

fi = a parameter in equation set 4-23 

H = liquid viscosity 

4, = heated perimeter 

h 

0^ = liquid density 

= vapor density 
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a 


surface tension 


0 = an experimentally evaluated bubble contact angle function 

DIFFERENTIALS 


dz 


two phase pressure gradient 


295 



APPENDIX C: ILLUSTRATIVE HEAT TRANSFER AND 
PRESSURE DROP DERIVATIONS 


1. Linear, Nonwetting Slug Flow Heat Transfer 

It is postulated that alternate slugs of liquid and vapor flow through 
the boiler tube; the vapor slugs have dry walls. It is presumed that 
heat transfer from the boiler tube to the liquid and vapor slugs is 
defined in terms of the wall-to-saturation temperature difference 
and liquid and vapor convective conductances. It is postulated that 
the mixed mean or bulk temperature of a vapor slug is equal to the 
saturation temperature. The mean heat transfer conductance is 
then defined as 


h(x) = h (x) y (x) + h (x) y (x) 
I 1 v v 


(a) 


where h^ and h^ are the liquid and vapor slug conductances, re- 
spectively, y f and y^ are the liquid and vapor volume fractions 
respectively, and x is the vapor quality. 

The relation between vapor quality and the volume fraction for this 
flow type is. 




PAf; E BUNK 


Wfy r nun® 
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X = 


1 + 


W V 

I y l 

w y 
v v 


1 



(b) 


where and are the liquid and vapor weight densities, respectively. 

The convective conductance for the liquid slugs in Equation (a) can 
be expressed in terms of the classical Nusselt viscous flow equation 
with uniform wall heat addition. 


fc _ 48 \ 

h f 11 D 


(c) 


where k is the thermal conductivity of the liquid and D is the tube 
x 

diameter. 


The convective conductance for a vapor slug under viscous flow is 


h 

v 


48 \r 

11 D 


(d) 


The corresponding function for turbulent flow is 


h 

v 



0.8 


= h v (x = l)(x)° 8 


<e) 
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where 


h (x 
v 


k 

= 1) = 0.023~~ Re 
D v 


0.8 


(x = 1 ) Pr y 


0.4 


(f) 


The quantity is the vapor flow rate, is the vapor thermal conductivity, 
Re y is the vapor Reynolds modulus and Pr^ is the vapor Prandtl modulus. 

The substitution of Equations (b) to (f) into Equation (a) yields 



for turbulent gas flow. 

The corresponding substitutions yield 



for viscous gas flow. 
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Equations (g) and (h) have been evaluated for properties that are 
representative of mercury in Figure 3-8 of Chapter 3 . 


2 . Annular Flow Pressure Drop (for Viscous -Turbulent Flow) 

Consider the case of a viscous annulus contiguous to the tube wall with 
an inner turbulent core. The pressure drop for the turbulent core 
( 0< r <r x > is, 


_ 

dx 




(i) 


where 


r 

r ] 

P 

x 

w 

V 


radial distance 

radius of the turbulent core 

pressure 

axial distance 

Weisbach friction factor 

density of turbulent fluid 

average velocity of turbulent fluid with respect to the moving 
interface 


The mass flow rate of the turbulent core is 


G 


1 



w x ( Vi + V) 


(j) 
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where the velocity of the moving Interface at r = r Is v^ . 


The quantity v is obtained by equating the shear stress-pressure drop 
function. 


T _ JLP. £ 
«x 2 

to the shear stress-shear strain function. 



(k) 


(I) 


The result is 


v 


1 



(m) 


where 

u , viscosity of viscous fluid 

o 

r , inside pipe radius 

o 
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The substitution of Equations (j) and (m) into Equation (i) yields 


He. 

dx 




dx, 


TP 



4r x g 




_L HE. 

n dx 
o 





trr. 


w 


1 


+ 




(n) 


If the turbulent fluid in the core fills the whole pipe, the pressure 
drop function would be 



where the subscript lp 


B_J 

2 

4 IT gr w 
B o 1 


refers to the inner turbulent fluid. 


(o) 


302 



Upon dividing Equation (n) by Equation (o) there results 




% 

where 


G = 
o 




2 2 9 

^ ( r o F 1 VI 

4» o dx \ 2 ~ 2 )\ 


(q) 
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Upon equating the pressure drop terms in die two fluids, a complicated 

algebraic equation for r /r results (not shown here) in terms of £ and 

1 o 

the dimensionless moduli 


gr /l 1 

q 


Gw G ii w 

o 1 , 1 1 o 

- , and 


G^w 


G #l w 
o o l 


thus, the two phase flow pressure drop ratio in Equation (p» is also 
a function of these moduli. 
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APPENDIX D: COMMENTS ON SOME ASSOCIATED TOPICS 


There are a number of associated technical topics that the boiler designer 
should be aware of, and several of them are briefly outlined below: 

1. Fluid Flow Stability in Boilers 

The boiler designer should have some understanding of fluid flow 
stability in single and multitube boilers. The flow dynamics in both 
cases are related to such parameters as nonuniform tube heating, 
liquid superheating, differences between the densities of the liquid 
and vapor, differences in pressure drop with flow regime type, and 
others . 

On a practical basis, flow restrictors in the form of orifices are used 
in multitube systems to assist in balancing the nonuniform flow patterns. 
There are a number of references in the literature for the reader if 
further information on this topic is desired (see references 1 to 4, 
for exam ile) . 

2. Nonequilibrium Effects 

There are some working fluids that can sustain significant amounts of 

liquid superheat; this phenomenon can occur even in forced flow systems. 

One working fluid that exhibits a significant amount of superheat prior 

5 

to flashing into steady state boiling is potassium. Although this problem 
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is generally not a critical one, the boiler designer should be cognizant 
of the fact that a sudden change in boiling state can occur with some 
working fluids, resulting in sudden changes in pressure and fluid flow. 
Likewise, techniques have been developed to overcome this phenomenon 
wher. it does occur such as the use of boiling initiators, artificial cavi- 
ties, and local high temperature surfaces. 

3. Working Fluid Stability 

The power plant designer must make a choice of the working fluids 
to be used in his system; the choice is usually based on a number of 
considerations. Ideally, it is desirable to have a working fluid that 
has good thermal properties so that optimum heat transfer performance 
can be obtained. In addition, it is necessary that the thermodynamic 
properties are also optimum (to yield good cycle performance) - For 
many applications, the working fluid should be nontoxic, nonflammable, 
noncorrosive, and of low cost. 

A final criterion is that the working fluid should have good chemical 
stability over the range of operating temperatures so that the fluid 
properties do not change with time or so that the decomposition products 
do not foul the prime heat transfer surfaces. It is suggested that the 
designer consult references on fluid decomposition, such as 6 and 7, 
for further information on this topic. 

4. "Boiler Conditioning" 

There are certain working fluids (primarily the liquid metals) which 
can exhibit heat transfer performance changes with time. Such fluids 


306 



must operate in very clean environments and the fluids themselves 
must be fairly pure- Oxygen, nitrogen, and lubricating oil leaks 
must not occur; if such foreign materials do enter the system, the 
surface physics and chemistry can change significantly so that a 
normally wetting fluid-wall system can be transformed to a partially 
or totally nonwetting system, thereby changing phase distributions 
and boiler performance. It is suggested that the boiler designer 
consult reference 8 which provides more information on this topic. 
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